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Preface 


This  report  summarizes  twelve  months  of  experimental  effort 
devoted  to  the  analysis  of  diatomic  molecules  consisting  of  one  Column 
1  atom  (alkali  metals)  and  one  Column  II  atom  (alkaline  earth).  This 
study  was  an  integral  part  of  the  research  efforts  of  the  New  Laser 
Concepts  Branch  of  the  Air  Force  Weapons  Laboratory,  Kirtland  AFB,  New 
Mexico.  The  Air  Force  Weapons  Laboratory  has  been  interested  in  these 
molecules  as  the  optically  active  medium  for  potential  laser  candidates. 
This  research  resulted  in  the  first  observation  and  analysis  of  the 
LiMg  and  LiCa  molecules. 

Several  different  types  of  experiments  were  performed  during  the 
course  of  this  work.  These  included  laser  induced  chemiluminescence 
experiments  for  observing  the  molecules  and  characterizing  their 
potential  energy  curves  and  fluorescence  time  delay  and  efficiency 
measurements  for  determining  reaction  rates.  Through  this 
experimentation  and  the  study  of  the  associated  theory,  I  have  been 
able  to  increase  my  understanding  and  appreciation  of  chemical  physics 
research  while  making  an  original  contribution  to  science. 

I  would  like  to  thank  my  advisor  Dr.  D.  J.  Benard  of  the  Air 
Force  Weapons  Laboratory  for  his  guidance  through  the  course  of  this 
project.  I  would  like  to  express  my  appreciation  to  Dr.  S.  J.  Davis 
for  his  advice  and  the  use  of  his  equipment  during  these  efforts  and 
Dr.  H.  Michels  of  United  Technology  Research  Center  for  his  timely 


ab  initio  calculations.  Additionally,  I  wish  to  thank  Dr.  A.  Hunter 
and  Dr.  D.  Shankland  and  the  remainder  of  my  committee  for  their 
participation  in  this  project.  I  would  like  to  acknowledge  Mr.  Gary 
Lee  for  his  help  in  performing  the  experiments.  Finally,  I  must 
acknowledge  the  three  people  without  whom  I  never  would  have  completed 
this  effort,  my  wife  Diane,  daughter  Nicole  and  son  Michael.  Without 
their  support,  encouragement  and  sacrifices  this  work  would  not  have 
been  possible. 


DAVID  K.  NEUMANN 
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Abstract 

Emission  from  LiCa  and  LiMg  excimers  was  obtained  for  the  first  time  upon 
pumping  mixtures  of  Li  and  Ca  vapors  and  Li  and  Mg  vapors  with  an  Ar+ 
laser.  The  emission  spectra  were  analyzed  spectroscopically  and  found  to 
be  in  good  agreement  with  a  synthetic  spectrum  derived  from  ab  initio 
calculations  of  the  LiMg  and  LiCa  potential  energy  curves. 

A  kinetic  model  was  developed  ta  explain  the  LiMg  formation  process. 

The  kinetics  of  production  of  LiMg*  were  analyzed  by  time  delay 
measurements  and  efficiency  measurements  of  the  fluorescence  intensity. 
The  LiMg*  emission  results  from  the  following  reactions: 

Li2  (X)  +  hu  (488  nm)  •>  Li2(B) 

Li2(B)  +  Mg  -*■  LiMg*  +  Li 
LiMg*  -*  LiMg  +  hu  (960  run) 

The  LiMg  formation  rate  constant  for  the  above  reaction  was  determined 
to  be  2.9  x  10  ^  cm3/sec.  Several  other  rate  constants  were  measured 
for  the  Li/Mg  system  which  resulted  in  an  overall  specification  of 
the  kinetic  model. 
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Chapter  I  -  Introduction 


A.  Excimers 

This  thesis  reports  the  first  study  of  the  kinetic  and  spectro¬ 
scopic  properties  of  gas  phase  LiMg  and  LiCa  diatomic  molecules. 

These  species,  composed  of  an  alkali  metal (Group  IA)  atom  and  an 
alkaline  earth(Group  IIA)  atom,  are  part  of  a  larger  class  of  diatomic 
molecules  known  as  excimers.  Excimers  are  characterized  by  an  un¬ 
stable  or  repulsive  molecular  ground  state  as  shown  in  Figure  1. 

Thus,  an  excimer  molecule  only  exists  as  a  stable  species  in  bound 
excited  states. 

The  noble  gas  halides  are  a  class  of  diatomic  molecules  observed 

1-4 

to  have  these  characteristics.  The  repulsive  nature  of  the  ground 

electronic  state  in  these  species  is  a  direct  result  of  the  closed 
shell  structure  of  the  noble  gas  atoms.  Excitation  of  the  noble  gas 
atoms  results  in  an  open  electronic  structure  with  unpaired  electrons. 
Consequently,  ionic  and  covalent  interactions  of  the  excited  noble  gas 
atoms  with  halogen  atoms  can  result  in  bound  molecular  states. 

Other  excimer  species  include  the  noble  gas  dimers,  noble  gas 
alkalis  and  the  alkaline  earth  dimers.'*  ^  In  the  latter  case,  the 
alkaline  earth  atoms  act  as  psuedo  noble  gas  species  since  their 
ground  electronic  states  consist  of  a  core  and  two  paired  s  electrons. 
Similar  considerations  apply  to  molecules  such  as  LiMg  which  are 
composed  of  an  alkali  metal  and  alkaline  earth  atom.  The  contribution 
to  bonding  due  to  ionic  interactions  is  small  because  the  ionic  states 
are  far  removed  in  energy.  There  is  little  potential  for  covalent 
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Figure  1  -  Potential  Energy  Curve  of  Excimer  Molecule 
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bonding  due  to  the  closed  shell  structure  of  the  alkaline  earth  atom. 
Therefore,  it  is  expected  that  LiMg  and  LiCa  in  particular  will  have 
the  repulsive  ground  states  that  are  characteristic  of  excimers. 

The  characteristic  emission  spectra  of  excimer  molecules  are 
strongly  affected  by  the  dissociative  nature  of  the  ground  state.  The 
bound  states  of  diatomic  molecules  are  characterized  by  discrete 
vibrational  and  rotational  energy  levels  resulting  from  the  quantized 
motion  of  the  nuclei.  In  an  excimer,  transitions  from  bound  excited 
states  terminate  in  a  continuum  of  energy  levels  characteristic  of  the 
repulsive  ground  state.  Consequently,  the  emission  spectrum  of  an 
excimer  is  continuous  as  is  shown  in  Figure  1. 

Detecting  the  existence  of  an  excimer  requires  the  use  of 
specialized  techniques  because  of  the  repulsive  ground  state. 

Standard  experimental  techniques  such  as  absorption  or  direct  optical 
excitation  will  work  only  under  conditions  where  a  large  portion  of 
the  ground  state  atoms  are  thermally  associated.  For  most  metal  vapor 
excimers  this  requires  very  high  densities  which  are  experimentally 
difficult  to  obtain  due  to  the  reactivity  of  the  alkali  metals. 
Additionally  in  high  density  regimes,  collisional  phenomena  such  as 
chemical  reaction,  electronic  quenching  and  energy  transfer  which 
compete  with  the  radiative  processes  may  inhibit  observation  of  the 
excimer.  In  the  experiments  conducted  for  this  thesis  a  laser  was 
used  to  optically  excite  an  alkali  metal  dimer  molecule  maintained  at 
low  densities.  Following  reaction  of  the  excited  dimers  with  alkaline 
earth  atoms  radiative  decay  of  the  excimer  was  observed.  This  process 
is  termed  laser  induced  chemiluminescence. 
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A  plot  of  the  reaction  symmetries  for  a  generic  alkali/  alkaline 
earth  system  is  shown  in  Fig.  2.  Analysis  suggests  that  laser  induced 
chemiluminescence  experiments  could  produce  observable  excimer  emission 
since  the  formation  of  the  excimers  via  the  reaction  of  excited  alkali 
dimers  with  alkaline  earth  atoms  is  both  symmetry  allowed  and  exothermic. 
Reaction  of  the  alkaline  earth  atom  with  the  ground  state  dimers  is  an 
endothermic  process.  However,  excitation  of  the  dimer  substantially 
reduces  its  binding  energy  and  provides  the  overall  energy  necessary 
to  make  the  excimer  production  reaction  energetically  feasible.  This 
demonstrates  that  different  states  of  the  same  chemical  species  can 
have  grossly  different  reactive  properties.  Thus,  each  state  must  be 
treated  as  a  separate  entity  in  a  kinetic  analysis.  Consequently,  the 
kinetics  of  even  this  relatively  simple  system  can  become  quite  complex. 
The  diagram  shows  that  the  excimers  can  also  be  produced  via  an  alternate 
mechanism  involving  the  reaction  of  ground  state  alkali  dimers  with 
metastable  alkaline  earth  atoms.  This  process,  however,  will  not  be 
examined  in  this  thesis. 

The  success  of  a  laser  induced  chemiluminescence  experiment  not 
only  depends  upon  the  existence  of  a  proper  kinetic  channel  for  the 
excimer  formation  but  also  upon  the  relative  rates  of  the  competing 
processes  that  consume  the  excited  alkali  dimer  molecules.  These 
include  rapid  radiative  decay  and  collisional  energy  transfer  to 
excited  atomic  states.  Additionally,  competition  between  the  radiative 
decay  of  the  excimer  molecule  and  chemical  reactions  of  the  excimer 
with  alkali  metal  atoms  can  substantially  limit  the  yield  of  excimer 
photons  obtained.  Consequently,  it  is  necessary  to  perform  the  laser 
induced  chemiluminescence  experiments  in  a  density  regime 
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that  is  favorable  to  both  the  chemical  formation  and  radiative  decay 
of  the  excimers.  This  requires  the  coincidence  of  a  high 
concentration  of  alkaline  earth  atoms  and  low  concentration  of  alkali 
metal  atoms.  These  requirements  tend  to  be  contradictory  since  the 
vapor  pressure  curves  for  the  alkali  metals  and  alkaline  earths  follow 
one  another.  A  compromise  must  be  made  to  achieve  optimal  excimer 
emission.  Fortunately,  Mg  and  Ca  are  more  volatile  than  Li  so  that  a 
reasonable  environment  for  LiMg  or  LiCa  laser  induced 
chemiluminescence  can  be  obtained  in  a  simple  heat  pipe. 

B.  Historical  Background 

The  first  observation  of  an  alkali  metal/alkaline  earth  diatomic 
molecule  occurred  long  before  the  concept  of  an  excimer  had  been 
postulated.  In  1898,  Liveing  and  Dewar  observed  diffuse  absorption 

g 

bands  in  high  pressure  cells  containing  Na  or  K  and  Mg.  In  1925, 
Barat  performed  essentially  the  same  experiments  in  which  he  recon¬ 
firmed  the  observations  of  NaMg  and  KMg  and  additionally  observed  RbMg 
9 

and  CsMg.  These  efforts  demonstrated  the  existence  of  gas  phase 
diatomic  species  composed  of  alkali  metal  and  alkaline  earth  atoms, 
but  provided  no  information  on  their  molecular  structure.  Following 
this  work,  there  were  no  efforts  concerning  IA/IIA  diatomics  until 
1976. 

In  the  mid  1970's,  research  efforts  investigating  gas  phase 
chemiluminescent  reactions  were  conducted,  with  the  ultimate  goal 
being  the  development  of  a  short  wavelength  chemical  laser.  Research 
in  this  area  by  Benard  and  Slafer  resulted  in  the  discovery  of  a 
chemical  chain  reaction  that  efficiently  produced  alkaline  earth  atoms 
in  a  metastable  electronically  excited  state. ^  Additional 
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experiments  were  performed  to  determine  if  the  energy  stored  in  the 

alkaline  earth  atoms  could  successfully  be  transferred  to  a  lasing 

species.  During  this  work,  Mg  metastables  were  reacted  with  dimers 

in  a  low  pressure  flow  tube.  Chemiluminescence  was  observed  which  was 

attributed  to  the  KMg  excimer.^  Subsequent  experiments  involving 

optical  excitation  of  in  a  static  cell  containing  low  pressures  of 

1 2 

K  and  Mg  again  resulted  in  observation  of  KMg  emission.  Later,  in 
unpublished  work  at  the  Air  Force  Weapons  Laboratory,  the  excimer  NaMg 
was  observed  in  an  analogous  manner.  A  variation  of  this  experiment 
in  which  K  was  added  to  the  Na  and  Mg  demonstrated  enhanced  NaMg 
emission.  Using  a  tunable  dye  laser  NaMg  emission  intensity  versus 
laser  wavelength  was  recorded.  The  intensity  profile  matched  the 
predicted  absorption  profile  for  the  NaK(X1I->C1n)  transition.  The 
result  established  that  the  excimer  emission  resulted  from  chemilumi¬ 
nescence  in  which  optically  excited  NaK(C^Fl)  was  the  reactive  pre¬ 
cursor.  Therefore,  the  emitting  species  was  not  directly  pumped  by 
the  laser.  This  experiment  originated  the  phenomena  of  laser  induced 
chemiluminescence . 

These  experiments  demonstrated  that  alkali  metal/alkaline  earth 
excimers  could  be  observed  in  emission  at  low  pressure.  However,  the 
results  of  these  experiments  were  qualitative  in  nature.  The  emitting 
species  and  precursor  had  been  identified,  but  no  detailed  investiga¬ 
tion  of  the  spectroscopy  or  kinetics  had  been  conducted. 

C.  Applications 

The  IA/IIA  excimers  observed  by  Benard  were  considered  as  chem¬ 
ically  pumped  excimer  laser  candidates.  The  potential  for  such  a 
laser  was  limited  in  the  case  of  NaMg  and  KMg  by  the  coincidence  of 
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the  excimer  emisssion  bands  and  their  respective  alkali  dimerCNa^  and 
I^)  absorption  bands.  Two  additional  factors  prevented  an  evaluation 
of  their  laser  potential:  kinetics  data  was  not  available  to  assess 
the  excimer  formation  or  loss  rates,  and  the  scalability  of  a  source 
of  metastable  alkaline  earth  atoms  had  not  been  demonstrated. 

If  a  new  alkali  metal/alkaline  earth  excimer  could  be  discovered, 
which  did  not  have  the  dimer  absorption  problem,  and  if  the  formation 
and  loss  rates  could  be  measured,  then  this  species  could  be 
considered  for  possible  use  in  a  laser  system.  At  the  very  least,  it 
might  have  potential  as  an  optically  excited  excimer  laser. 

D.  Research  Objectives 

The  research  objectives  for  this  thesis  were: 

(a)  to  investigate  a  new  alkali  metal/alkaline  earth 
excimer  using  the  technique  of  laser  induced  chemiluminescence 

(b)  to  characterize  the  molecular  properties  of  the 
observed  emitter  using  available  experimental  and  theoretical  infor¬ 
mation 

(c)  to  determine  the  kinetic  processes  and  rates  important 
in  the  production  and  decay  of  the  excimer. 

E.  Accomplishments 

Initially,  various  alkali  metal,  alkaline  earth  combinations  were 
evaluated  using  the  laser  induced  chemiluminescence  technique.  The 
combinations  of  Li  and  Mg  and  Li  and  Ca  resulted  in  the  observation  of 
new  emissions  in  the  900-1100  nm  region  which  were  attributed  to  the 
respective  excimers,  LiMg  and  LiCa.  The  emission  spectra  of  the  mole¬ 
cules  were  recorded,  analyzed  and  compared  to  the  results  of  ab  initio 
calculations  completed  by  H.  Michels  at  the  United  Technology  Research 


Center.  Using  the  results  of  the  emission  experiments  and  theoretical 
considerations,  a  kinetic  model  for  the  LiMg  system  was  postulated. 
Experiments  were  then  conducted  to  investigate  the  kinetics  of  the 
LiMg  system.  Modulation  was  imposed  on  the  kinetic  system  by  electro- 
optically  chopping  the  laser  excitation  source.  Time  delays  be¬ 
tween  modulated  signals  from  various  emitters  were  measured  and  re¬ 
lated  to  kinetic  rates.  Additionally,  the  efficiency  of  light  produc¬ 
tion  was  measured.  The  results  of  these  experiments  were  then  used  to 
quantitatively  determine  kinetic  rates  for  the  postulated  model. 

A  total  of  ten  reaction  rate  constants  and  three  radiative  lifetimes 
were  calculated,  measured  or  had  upper  bounds  established  as  a 
result  of  the  kinetics  work.  The  kinetics  work  together  with  the 
analysis  of  emission  from  two  new  excimers  constituted  the  thesis 
work. 
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Chapter  II  -  Spectroscopy  of  LiCa  and  LiMg 


A.  Introduction 

This  thesis  reports  the  first  observation  of  the  diatomic  species 
LiMg  and  LiCa.  Although  absorptions  in  several  alkali-Mg  diatomics 
were  seen  before  no  alkali-Ca  diatomic  had  ever  been  observed.  The 
work  presented  in  this  chapter  represents  the  first  attempt  to  characterize 
LiMg  and  LiCa  excimers.  Experimental  evidence  of  the  existence  of 
these  species  is  provided.  A  qualitative  analysis  of  the  electronic 
structure  of  LiMg  and  LiCa  is  presented.  The  results  of  molecular  ab 
initio  calculations  including  the  potential  energy  curves  and  predicted 
emission  spectra  for  LiMg  and  LiCa  are  explained  and  related  to  the 
experimental  work.  The  theoretical  calculations  are  the  unpublished 
work  of  Dr  H.  Michels  of  UTRC. 

B.  Theory 

1.  Molecular  State  Characterization 

The  characteristics  of  the  ground  state  and  first  few  excited 
states  of  LiMg  and  LiCa  can  be  qualitatively  predicted  from  knowledge 
of  the  electronic  structure  of  the  separated  atoms.  A  ground  state  Li 
atom  (2S)  has  a  single  s  electron  in  its  outer  valence  shell  while  the 
ground  state  alkaline  earth  atoms  (1S)  have  two  paired  s  electrons. 

The  first  excited  state  of  the  Li  atom  (2P)  results  from  the  promotion 
of  the  s  electron  to  a  p  orbital  (p^,  p^,  Pz)-  The  first  excited 
state  of  the  alkaline  earth  atoms  (3P)  also  results  from  promotion  of 
an  s  electron  to  a  p  orbital.  The  Li(2P)  state  is  connected  to  the 
ground  state  by  a  highly  allowed  electric  dipole  transition.  The 
Ca(3P)  and  Mg(3P)  states  lie  higher  in  energy  than  the  Li(2P) 
state  and  transitions  to  their  respective  ground  states  are 
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spin  forbidden.  An  energy  level  diagram  for  Li,  Ca ,  and  Mg  is  shown  in 

Figure  3a.  The  molecular  states  correlating  to  these  separated  atom 

1 3  *“  1 6 

limits  have  been  derived  from  symmetry  principles  and  are  shown 

for  LiMg  in  Figure  3b.  Since  the  energy  ordering  of  the  separated 
atom  limits  are  the  same  for  both  LiMg  and  LiCa,  the  following  discus¬ 
sion  applies  to  LiCa  as  well. 


20 
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Figure  3a  -  Li,  Ca ,  and  Mg  Atomic  Energy  Levels 
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Figure  3b  -  LiMg  Correlation  Diagram 

LiMg  has  a  single  2I  ground  electronic  state  correlating  to  the 
Li(2S)  and  Mg(xS)  separated  atom  limits.  The  first  manifold  of 
excited  states  are  of  2I  and  2n  symmetry  and  correlate  to  ground  state 
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Mg(*S)  and  excited  Li(2P)  atoms.  A  manifold  of  four  electronic  states 
of  symmetry  41 ,  4n,  2I  and  2I1  correlate  to  excited  Mg(3P)  and  ground 
state  Li(2S)  atoms. 

Some  aspects  of  the  molecular  bonding  in  these  electronic  states 

can  be  predicted  from  consideration  of  the  atomic  configurations  and 

the  resulting  interatomic  forces.  The  2X  ground  state  is  expected  to 

be  predominantly  dissociative  in  character  due  to  short  range 

repulsive  valence  forces.  These  forces  which  consist  of  both 

coulombic(charge)  and  exchange (Paul i  Principle)  contributions,  result 

from  the  interaction  of  the  outer  s  electron  of  the  Li  atom  with  the 

paired  s  electrons  of  the  Mg  atom.  At  longer  internuclear  distances, 

it  is  expected  that  the  repulsive  forces  will  be  dominated  by  weak 

attractive  dispersion  forces  resulting  in  a  small  potential  minimum. 

Dispersion  forces  result  from  interactions  of  induced  multipole 

moments  of  the  atoms  and  usually  are  dominated  by  the  induced  dipole- 

induced  dipole  interaction  which  is  more  commonly  called  the  Van  der 

Waals  force. ^  Since  the  induced  dipole  moment  is  directly 

proportional  to  the  polarizability  of  the  atoms,  it  is  expected  that 

1 8 

the  larger  polarizability  of  Ca  will  result  in  a  more  strongly  bound 
ground  state  for  LiCa  than  for  LiMg. 

19  42 

As  in  the  case  of  LiHe  and  NaHe  ’  the  bonding  in  the  first 
set  of  excited  states  can  be  directly  related  to  the  three  possible 
orientationsCp^  or  p^,  p^)  of  the  excited  Li  p  orbital.  With 
the  electron  in  a  pz  orbital,  aligned  to  the  internuclear  axis  of  the 
diatomic  molecule  and  interacting  with  the  paired  outer  s  electrons  of 
the  alkaline  earth  atom,  the  repulsive  forces  should  be  even  stronger 
than  for  the  ground  state  LiCa  and  LiMg  interactions.  The  resulting  2X 


state  should  be  strongly  repulsive.  The  2Ii  state  should  be  bound  due 
to  attractive  dispersion  forces  (ion-induced  dipole)  between  the  un¬ 
screened  Li+  core  resulting  from  the  p^  or  p^  configuration  of  the 
outer  electron  and  the  Mg  outer  s  electrons. 

The  set  of  excited  2I  and  2n  states  correlating  to  Mg(3P)  and 
Li(2S)  should  be  strongly  bound  due  to  attractive  valence  forces 
resulting  from  interactions  of  the  unpaired  Mg  electrons  with  the 
single  Li  s  electron.  As  discussed  for  the  previous  set  of  2I  and  2n 
states,  the  2fl  state  should  be  more  strongly  bound  due  to  the 
additional  attractive  ion-induced  dipole  dispersion  forces. 

The  quartet  states  in  this  manifold  should  be  strongly  repulsive 
due  to  large  repulsive  exchange  forces  resulting  from  the  interaction 
of  spin  aligned  electrons. 

The  bonding  exhibited  by  the  excited  electronic  states  of  LiMg 

and  LiCa  is  not  only  affected  by  forces  resulting  from  the  atom-atom 

interactions,  but  by  the  perturbation  of  one  state  by  another  of  like 

symmetry.  The  upper  set  of  LiMg  and  LiCa  2fl  and  21  states  are 

expected  to  be  strongly  bound  while  the  lower  set  is  primarily 

repulsive.  This  normally  would  result  in  crossings  of  the  potential 

energy  curves.  However,  the  potential  energy  curves  of  states  of  the 

same  symmetry  cannot  cross  in  the  adiabatic  limit  but  tend  to  repel  one 
20  21 

another.  ’  The  net  effect  for  LiMg  and  LiCa  should  be  strongly 

bound  lower  ZT1  and  2I  states  and  repulsive  upper  2n  and  2I.  The 
effect  of  the  interaction  should  be  more  pronounced  in  LiCa  because  of 
the  closer  resonance  of  the  Li(2P)  and  Ca(3P)  states  which  should 
result  in  stronger  mixing  of  the  wavefunctions . 
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Based  on  this  analysis,  these  states  can  now  be  given  the  standard 
spectroscopic  notation  according  to  their  energy  ordering.  The  ground 
state  will  be  labeled  X21,  the  first  two  excited  states  will  be  labeled 
A2n  and  B2!  and  the  next  two  states  C21I  and  02I. 

The  repulsion  of  the  2fl  and  21  states  should  significantly  affect 
the  observed  transitions  in  LiMg  and  LiCa.  The  lower  set  of  states 
are  stable  enough  to  support  substantial  concentrations  of  molecules 
given  an  appropriate  chemical  production  mechanism.  The  upper  set  of 
states  should  be  repulsive  and  therefore  will  not  support  any  concentrations. 
Observed  emissions  in  LiMg  and  LiCa  should  therefore  result  from 
allowed  transitions  between  the  lower  2fl  and  2Z  states  and  the  ground 
state . 

A  quantitative  determination  of  the  electronic  structure  of  LiMg 

and  LiCa  has  been  obtained  by  means  of  ab  initio  calculations  of  the 

24 

potential  energy  curves.  These  calculations  were  completed  by  H. 

Michels  of  United  Technology  Research  Center  as  part  of  a  cooperative 

effort  to  obtain  both  a  theoretical  and  experimental  understanding  of 

these  excimers.  The  calculational  procedure  that  was  utilized  was 

22 

the  valence  configuration  interaction  (VCI)  method.  For  these 

calculations  a  spin  free  nonrelativistic  Hamiltonian  was  employed  in 

the  Born-Oppenheimer  approximation.  Symmetry  compatible  electronic 

wavef unctions ,  composed  of  linear  combinations  of  atomic  orbitals  and 

depending  on  the  nuclear  coordinates  were  made  to  be  optimum  approximations 

to  solutions  of  the  Schrodinger  equation  by  involving  the  variational 
23 

principle . 

A  full  VCI  expansion  was  chosen  for  the  LiCa  studies  using  an 
optimized  minimum  basis:  Li(ls),  Ca(ls,  2s,  2p)  for  the  core 
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electrons  and  a  double  zeta  optimized  basis:  Li(2s,2s,  2p,2p), 

Ca(3s,  3s,  3p,  3p)  for  the  valence  shell  description.  This  basis 
yielded  a  2S-2P  Li  excitation  energy  of  14830/cm  (14904/cm  =  experi¬ 
mental  value)  and  a  1S-3P  Ca  excitation  energy  of  15135/cm  (15225/cm  = 
experimental  value),  thus  insuring  a  proper  long-range  mixing  of  the 

ry  / 

wavefunctions  connecting  to  Li*(2P)+Ca ( *S)  and  Li (2S)+Ca*(3P) -  A 

similar  basis  set  was  used  for  the  LiMg  studies  with  similar  results. 

27 

The  calculated  potential  energy  curves  for  LiCa  and  LiMg  are 
shown  in  Figures  4  and  5.  As  predicted,  qualitatively  the  ground 
states  of  both  molecules  exhibit  small  van  der  Waals  potentials.  The 
LiCa  ground  state  is  more  strongly  bound  (580/cm)  than  the  LiMg  ground 
state  (137/cm).  The  A2F1  and  B 2I  states  are  more  strongly  bound  than 
the  C  and  D  states.  The  coefficients  of  the  wavefunction  expansions 
showed  that  the  excited  electronic  states  of  both  molecules  were 
strongly  mixed  even  at  relatively  large  internuclear  separation.  At 
shorter  internuclear  separation,  the  A2n  and  BzZ  states  assumed  more 
of  the  character  of  the  alkaline  earth  (3P)  and  Li(2P)  states 
respectively.  The  changing  characters  of  the  electronic  states 
with  internuclear  separation  strongly  affects  their  transition 

2  2 

probabilities  to  the  ground  state  because  of  the  different  Li(  P-  S) 

3  1 

and  Ca(  P-  S)  transition  probabilities.  The  magnitude  of  this  effect 

will  be  determined  in  the  next  section.  Calculations  of  the  quartet 

states  were  not  performed  since  they  do  not  optically  connect  to  the 

ground  state  and  do  not  interact  with  the  doublet  states. 

The  spectroscopic  constants  for  the  ground  state,  A  and  B 
27 

states  of  LiMg  and  LiCa  are  shown  in  Tables  1  and  2.  The  energy  levels 
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Table  I  -  LiCa  Calculated  Spectroscopic  Constants 


Spectroscopic 


State 


Constant 


X2I 


A2n 


B2! 


Te(ev) 

0.0000 

1.0800 

1.1700 

w^Ccm  *) 

93.3600 

252 .4400 

197.4400 

^XgCcm'1) 

7.3200 

1.2900 

.9600 

ae(cm  *) 

.0119 

.0019 

.0009 

B  (cm  *) 
e 

.1890 

.2740 

.1950 

r  (A) 

3.8830 

3.2260 

3.8220 

D  (eV) 
e 


.0720 


.8820 


.7650 


Table  II  -  LiMg  Calculated  Spectroscopic  Constants 


Spectroscopic 

State 

Constant 

X2I 

Azn 

B2! 

Te(ev) 

0.0000 

1.1200 

1.5619 

u)e(cm  *) 

65.7000 

302.5000 

328.2000 

u>exe(cm‘l) 

6.9000 

2.0000 

8.7000 

0e(cm'1) 

.0160 

.0032 

.0040 

Be(cm-1) 

.2024 

.4010 

.2770 

re(A) 

3.9280 

2.7900 

3.3560 

De(eV) 

.0167 

.7630 

.3000 

D0(eV) 

.0073 

.7450 

.2800 
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of  each  state  are  described  by  the  Dunham  expression: 

E(e,v,J)  =  Tg  +  uie(v+5i)  +  weXe(v+^)2  +'“ 

BgJ (J+i )  -  «e(v+y  (j(j+d) 

where  Tg  is  the  electronic  energy,  u)e  is  the  vibrational  spacing,  ^eXe 

is  the  anharmonic  correction,  B  is  the  rotational  constant  and  «  is 

e  e 

the  vibrational  and  rotational  coupling  constant.  The  remaining  terms 

in  Tables  1  and  2  are  rg  the  internuclear  separation  where  the 

potential  energy  curve  reaches  a  minimum  and  Dg  and  Dq  the 

dissociation  energies  of  the  state  measured  from  the  bottom  of  the 

potential  well  and  first  vibrational  level  respectively. 

2.  Electronic  Transitions  in  LiMg  and  LiCa 

Radiative  transitions  between  the  excited  2I1  and  2I  states  and 

the  2I  ground  state  of  LiMg  and  LiCa  are  electric  dipole  allowed 
29 

transitions.  The  transitions  in  LiMg  and  LiCa  are  expected  to  be 

both  bound-bound  and  bound-continuum,  based  on  application  of  the 

Franck-Condon  principle  to  Figures  4  and  5.  Electric  dipole  transi- 

30 

tions  in  bound- continuum  systems  have  been  analyzed  by  Mies  and  have 

31 

been  calculated  for  LiMg  and  LiCa  by  Michels.  The  bound  and 
continuum  spectra  are  characterized  by  matrix  elements  for  the  tran¬ 
sitions  where  the  vibrational  wavefunctions  are  obtained  by  numerical 

32  33 

solution  of  the  Schrodinger  equation  using  the  Numerov  method  ’  and 

the  continuum  wavefunctions  are  obtained  using  the  method  of  Allison 
34 

and  Dalgarno.  In  both  cases  the  dipole  moment  function  which  varies 
with  internuclear  distance  is  required  to  calculate  the  matrix 
elements.  The  dipole  moment  function  is  calculated  from  the  molecular 
wavefunctions  and  dipole  moment  operator  for  the  transition.  These 
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calculated  dipole  moment  functions  for  the  A2H,  B2!,  C2F1  and  D2!  to 
X2I  transitions  for  both  LiMg  and  LiCa  are  shown  in  Figures  6  and  7. 

At  shorter  internuclear  distances  the  21-21  transitions  have  stronger 
transition  moments  because  the  excited  state  2I  molecular  wavefunction 
is  dominated  by  the  Li(2P)  atomic  orbital.  The  Li(2P-2S)  transition 
has  a  large  transition  probability  (A=108/sec).  The  2n  states  are 
dominated  by  the  character  of  3P  states  of  the  alkaline  earth  atoms 
resulting  in  a  smaller  dipole  moment  function.  The  alkaline  earth 
(3P-1S)  transition  probabilities  are  small  (A=104/sec).  At  large 
internuclear  separation,  the  dipole  moment  functions  are  more  strongly 
affected  by  the  wavefunctions  associated  with  the  separated  atom 
limits.  The  A2n  and  B2!  states  have  larger  transition  moments  which 
result  from  the  dominance  of  the  Li(2P)  character  of  the  wavefunction 
while  the  C2n  and  D2!  tend  to  the  alkaline  earth  3P  limit  with 
resulting  smaller  transition  probabilities.  In  all  cases  the  2n 
states  retain  more  of  the  3P  character  resulting  in  smaller  overall 
dipole  moment  functions  for  the  2n-2I  transitions  than  for  the  corres¬ 
ponding  2Z-21  transitions. 

The  dipole  moment  functions  shown  in  Figures  6  and  7  were  used  to 

evaluate  bound-continuum  and  bound-bound  matrix  elements  and  to 

predict  the  emission  spectrum  resulting  from  A2FI-X2Z  and  B2£-X2I 
27 

transitions.  The  results  are  shown  in  Tables  3  through  10.  The 
bound-contiuum  emission  tables  for  LiCa(III,IV)  and  for  LiMg(VII ,VIII) 
show  the  center  wavelengths  of  various  peaks  in  the  emission  spectrum 
resulting  from  transitions  from  upper  state  vibrational  levels  to  the 
continuum  of  the  ground  state.  Since  each  level(v)  has  a  wavefunction 
with  v+1  nodes,  the  reflection  of  these  nodes  into  the  continuum 
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Figure  7  -  Dipole  Transition  Moments  for  LiCa 
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result  in  the  multiple  peaks  shown.  Not  all  of  the  peaks  for  all  the 
transitions  are  shown  since  the  transition  moments  become  negligible 
with  increasing  wavelength. 

The  bound-continuum  tables  also  list  the  relative  amplitudes  of 
the  transitions.  In  both  LiMg  and  LiCa  the  B-X  transitions  have  the 
largest  vibrational  transition  amplitudes.  This  is  due  to  the  larger 
transition  moments  for  the  1 -1  transitions  and  the  fact  that  the  B 
state  rg  are  larger  than  the  A  state  resulting  in  transitions  to 
flatter  portions  of  the  continuum  curve.  For  the  case  of  LiCa,  it  is 
apparent  that  most  of  the  bound-free  emission  should  be  localized  in 
the  region  of  1010-1020  nm  and  should  result  from  v'  =  3, A, 5  of  the  B 
state.  In  LiMg  the  major  contributions  to  the  bound-free  spectrum 
should  originate  from  the  B  state;  v*  =  0,1, 2, 3, 4, 5,  and  should  appear 
in  the  region  of  720-820nm. 

The  bound-bound  emission  tables  for  LiCa(V,VI)  and  LiMg(IX,X) 
show  that  the  strongest  emissions  are  the  1-1  emission.  Again  this 
can  be  attributed  to  both  the  stronger  transition  moments  and  the 
location  of  the  potential  wells.  In  LiCa  however  higher  vibrational 
levels  of  the  A  state  (v'=3,etc)  have  nonnegligible  transition 
amplitudes.  Overall  the  largest  transition  amplitudes  in  LiCa  result 
from  v'=3  in  the  A2n  states  and  v'=0,  1  in  the  B 2I  state.  In  LiMg  the 
major  contributions  to  the  bound  bound  spectrum  should  originate  from 
the  B2!  state  and  should  appear  in  the  region  of  690-750. 

The  emission  spectrum  resulting  from  the  C2n  and  D2Z  states  of 
LiMg  and  LiCa  to  the  X2I  ground  state  were  not  discussed  because  these 
states  are  primarily  repulsive  and,  therefore,  would  support  only  a 
very  small  steady  state  density  of  molecules.  Additionally,  these 


Table  III  -  LiCa (A2n-X2I)  Bound-Free  Emis  sion  Spectrum 


Upper  State 

Center 

Relative 

Vibrational  Level 

Wavelength(nm) 

Amplitude 

v'=0 

1237.5 

.010 

1 

1178.5 

.012 

1343.2 

.003 

2 

1145.0 

.001 

1205.0 

.007 

1394.9 

.002 

3 

1113.7 

.006 

1133.0 

.012 

1234.5 

.004 

1434.5 

.001 

4 

1084.5 

.017 

1154.0 

.008 

1266.0 

.003 

1486.0 

.000 

5 

1056.9 

.005 

1093.3 

.013 

1170.7 

.006 

1286.1 

.002 

6 

1047.7 

.021 

1103.7 

.010 

1182.6 

.005 

1300.5 

.002 
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Table  IV  -  LiCa (B2I-X2I)  Bound-Free  Emis  sion  Spectrum 


Upper  State 

Center 

Relative 

Vibrational  Level 

Wavelength  (run) 

Amplitude 

v '  —  0 

1091.4 

.000 

1 

1069.1 

.017 

1106.3 

.000 

2 

1055.9 

.045 

1102.8 

.001 

1154.2 

.000 

3 

1027.6 

.351 

1052.9 

.059 

1099.6 

.001 

1150.6 

.000 

4 

1008.2 

.907 

1012.3 

1.00 

1050.1 

.075 

1096.5 

.075 

1152.6 

.002 

5 

1009.0 

.907 

1047.4 

.088 

1093.6 

.003 

1149.3 

.000 

1211.1 

.000 
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Table  V  -  LiCa  (B2I-X2I)  Bound-Bound  Emiss  ion  Spectrum 
\(nm) ;  Relative  Intensity(Energy/sec) 


V 

v" 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1047.2 

1.00 

1026.3 
.  154 

1006.4 
.  172 

987.5 

.071 

969.4 

.049 

952.2 

.025 

935.7 

.014 

920.0 

.008 

1 

1053.2 

.529 

1032.1 

.075 

1012.0 

.035 

992.8 
.  140 

974.6 
.  122 

957.2 
.  128 

940.5 
.  101 

924.6 

.079 

2 

1059.3 

.183 

1037.9 

.452 

1017.6 

.050 

998.2 

.078 

979.8 

.002 

962.2 

.003 

945.4 

.029 

929.3 

.055 

3 

1065.7 

.029 

1044.1 

.589 

1023.5 

.013 

1004.0 

.002 

985.3 

.091 

967.5 

.068 

950.5 

.058 

934.3 

.019 

Table  VI  -  LiCa (A2FI-X2I)  Bound-Bound  Emission  Spectrum 
A(nm);  Relative  Intensity (Energy/sec) 


V 

v" 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1139.0 

.005 

1107.4 

.035 

1077.9 

.124 

1050.1 

.295 

1024.0 

.527 

999.5 

.761 

976.3 

.932 

954.4 
1 .00 

1 

1149.0 

.011 

1116.8 
.  064 

1086.8 
.  184 

1058.6 

.329 

1032.1 

.404 

1007.2 
.  346 

983.6 
.  189 

961.4 

.041 

2 

1157.7 

.018 

1125.1 

.095 

1094.6 

.219 

1066.0 

.289 

1039.2 

.220 

1013.9 

.072 

990.0 

.000 

967.6 

.095 

3 

1166.2 

.029 

1133.1 
.  126 

1102.2 

.229 

1073.2 

.207 

1046.0 

.070 

1020.4 

.001 

996.3 

.110 

973.5 

.288 
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Table  VII  -  LiMgCA2Fl-X2I)  Bound-Free  Emission  Spectrum 


Upper  State 

Center 

Relative 

Vibrational  Level 

Wavelength(nm) 

Amplitude 

v 1  —  0 

1338.2 

.006 

1 

1211.3 

.010 

1505.5 

.002 

2 

1127.3 

.016 

1291.7 

.004 

1610.6 

.010 

3 

1082.5 

.024 

1186.1 

.007 

1354. 6 

.002 

1686.2 

.001 

4 

1024.6 

.034 

1107.0 

.010 

1215.5 

.004 

1393.0 

.002 

5 

988.7 

.047 

1047.1 

.016 

1133.3 

.007 

1247.3 

.003 

6 

955.7 

.062 

1002.1 

.022 

1062.2 

.011 

1151.0 

.006 

7 

925.4 

.075 

968.8 

.031 

1016.4 

.016 

1159.0 

.009 
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Table  VIII  -  LiMg(B2I-X2n)  Bound-Free  Emiss  ion  Spectrum 


Upper  State 

Center 

Relative 

Vibrational  Level 

Wavelength(nm) 

Amplitudi 

< 

II 

o 

802.2 

.828 

1 

783.4 

.966 

841.8 

.552 

2 

757.3 

1.00 

771.4 

.448 

811.8 

.310 

3 

745.0 

.638 

754.0 

.552 

787.5 

.310 

813.3 

241 

4 

726.1 

.706 

748.0 

.310 

766.5 

.276 

786.0 

.241 

811.7 

.207 

5 

714.5 

.793 

763.8 

.259 

783.0 

.224 

808.5 

.207 

847.2 

.397 
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.  ible  IX  -  LiMg(A2n-X2I)  Bound-Bound  Emission  Spectrum 
A.(nm) ;  Relative  Intensity (Energy/sec) 


v' 

0 

1 

2 

3 

4 

0 

1111.0 

.000 

1075.4 

.000 

1042.5 

.000 

1011.9 

.000 

983.6 

.000 

1 

1117.1 

.000 

1081.2 

.000 

1047.9 

.000 

1017.0 

.000 

988.4 

.001 

2 

1121.9 

.000 

1085.6 

.000 

1052.0 

.000 

1020.9 

.000 

992.1 

.001 

v' 

5 

6 

7 

8 

9 

0 

957.2 

.001 

932.6 

.004 

909.7 

.013 

888.3 

.032 

868.1 

.072 

1 

961.7 

.003 

936.9 

.009 

913.9 

.023 

892.3 

.052 

871.9 

.099 

2 

965.2 

.003 

940.7 

.008 

917.0 

.020 

895.2 

.040 

874.8 

.068 

Table 

X  -  LiMg(B2I-X2I)  Bound-Bound  Emission  Spectrum 
A.(nm) ;  Relative  Intensity  (Energy/sec) 

v’ 

0 

1 

2 

3 

4 

0 

795.8 

.004 

780.5 

.024 

766.3 

.088 

753.2 

.228 

741.3 

.443 

1 

799.0 

.008 

783.5 

.042 

769.2 

.124 

756.1 

.240 

744.0 

.306 

v’ 

5 

6 

7 

8 

9 

0 

730.3 

.669 

720.3 

.807 

711.2 

.791 

702.9 

.634 

695.6 

.408 

1 

732.9 

.235 

722.9 

.069 

713.7 

.007 

705.4 

.224 

698.0 

.665 
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states  are  not  accessible  in  the  experiments  conducted  in  this  thesis 
due  to  energy  considerations  which  will  be  explained  in  the  discussion 
of  the  experimental  results. 

Finally  there  was  no  attempt  to  incorporate  the  effects  of 
temperature  and  state  dependent  chemical  production  rates  into  the 
calculations.  The  calculations  show  the  transition  amplitudes  derived 
from  quantum  mechanical  considerations.  To  relate  these  transition 
amplitudes  to  real  intensities  it  is  necessary  to  consider  how  an 
ensemble  of  molecules  is  distributed  among  all  the  possible  energy 
levels  as  a  result  of  collisional  effects.  These  considerations  will 
be  included  in  the  comparison  of  the  experimental  data  to  the  theory. 
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C.  Experimental  Considerations. 
1.  Heatpipe  Operation 


The  spectroscopic  experiments  performed  for  recording  emission 
spectra  from  LiMg  and  LiCa  required  metal  atom  densities  of  approxi¬ 
mately  1016/cm3.  This  required  heating  the  metals  to  temperatures  of 

600-1000°C.  In  1969,  Vidal  and  Cooper  developed  a  device  called  a 

36 

heatpipe  oven  for  use  in  spectroscopic  experiments,  which  had  the 
capability  to  maintain  uniform  densities  of  reactive  metal  vapors  over 
long  periods  of  time.  Figure  8  shows  the  design  of  the  heatpipe,  which 
consists  of  a  straight  pipe  sealed  at  both  ends  with  windows 
containing  ports  for  access  to  a  vacuum  pump  and  inert  gas  supply. 
Cooling  chambers  are  located  at  each  end  and  a  heater  coil  is  wrapped 
around  the  central  region  to  heat  the  metal  contained  inside.  Inside 
the  heated  region,  internal  to  the  pipe,  is  a  wire  mesh  material  which 
serves  to  reflux  condensed  metal  vapor  back  to  the  hot  interior.  When 
the  device  is  operated,  it  is  first  filled  with  an  inert  gas  and  the 
metal  to  be  vaporized.  When  heated,  the  metal  melts,  wets  the  wire 
mesh  wick  and  then  evaporates.  The  vapor  diffuses  toward  both  ends 
until  it  reaches  a  cooler  region  where  it  condenses.  The  condensed 
metal  flows  toward  the  center  heated  portion  of  the  pipe  as  a  result 
of  capillary  action  on  the  wick.  When  sufficient  heat  is  added  to 
achieve  a  temperature  such  that  the  metal  vapor  pressure  equals  the 
pressure  of  the  inert  gas  previously  introduced  into  the  cell,  the 
cycling  of  the  metal  sets  up  a  region  at  the  center  of  the  heatpipe 
where  only  the  metal  vapor  exists  at  a  vapor  pressure  equal  to  the 


inert  gas  pressure. 


Well  defined  regions  or  zone  boundaries  are  observed  between  the  metal 

vapor  and  the  inert  gas.  With  Li  in  the  heatpipe,  experiments  have 

demonstrated  that  the  temperature  is  uniform  within  the  heated  region 
36 

to  within  1°C-  When  two  metals  having  different  vapor  pressure 

curves  are  introduced  into  the  heatpipe,  a  second  set  of  zone 

boundaries  are  formed.  The  central  region  of  the  heatpipe  contains  a 

mixture  of  both  metals,  with  the  higher  vapor  pressure  metal  confining 

the  lower  vapor  pressure  metal.  The  second  region  contains  only  the 

37 

higher  vapor  pressure  metal  which  is  confined  by  the  inert  gas. 

When  both  metals  used  are  Group  I  metals  which  readily  wet  the 

wick  and  flow  back  to  the  center  of  the  heatpipe,  a  situation  such  as 

that  described  above  can  be  maintained.  Wicking  of  Group  II  metals  is 

much  more  difficult  to  obtain  and  in  many  cases,  the  heatpipe  must  be 

operated  with  the  buffer  gas  pressure  in  excess  of  the  metal  pressure 

to  better  confine  the  metal.  In  this  case,  no  well  defined  zone 

boundaries  are  formed  and  no  pure  metal  regions  are  formed.  The  inert 

gas  pressure  in  the  metal  region  is  just  the  difference  between  the 

inert  gas  pressure  measured  in  the  cold  region  and  the  metal 

pressure  defined  by  the  temperature.  Even  with  these  changes,  the 

same  processes  of  evaporation,  condensation,  and  flow  continue. 

Another  physical  parameter  that  must  be  considered  in  the 

operation  of  a  two  metal  heatpipe  is  the  solubility  of  one  metal  in 

another.  If  one  metal  is  completely  soluble  in  another,  the  vapor 

pressure  of  the  solution  is  reduced  by  the  molar  ratio  of  the  two 

38 

metals  (Raoult's  Law).  This  result  holds  for  dilute  two  component 
solutions.  Deviation  from  the  dilute  solution  situation  requires 
experimental  verification  tha*  metal  vapor  pressures  have  not  been 
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perturbed.  However,  the  lower  the  solubility  of  the  two  metals  in 

each  other,  the  smaller  should  be  the  effect  on  the  vapor  pressures. 

The  solubility  of  Group  II  metals  in  Group  1  metals  is  known  to  be  2% 
39-41 

or  less  which  should  result  in  small  deviations  to  the  vapor 

pressure  curves. 

2.  Apparatus 

The  mixed  metal  heatpipe  used  in  the  experiments  for  this 
thesis  was  modified  from  the  design  of  Cooper  and  Vidal  to  accommodate 
optical  excitation  of  the  metals  and  observation  of  the  resulting 
emission.  The  heatpipe  was  formed  in  the  shape  of  a  cross  as  shown  in 
Figures  9  and  10.  This  allowed  a  laser  beam  to  be  passed  through  in 
one  direction  and  the  resulting  fluorescence  to  be  viewed  from  the 
side.  The  design  was  slightly  modified  from  one  developed  by  D.  J. 
Benard  to  achieve  1000°C  operation.  The  modification  involved  removal 
of  all  brass  and  copper  parts  since  their  melting  points  were  being 
exceeded.  The  heatpipe  was  fashioned  from  a  solid  stainless  steel 
cylinder,  4"  in  diameter  and  6"  high.  The  cylinder  was  bored  to  form 
a  central  chamber  and  to  accept  four  lateral  sidearms  and  one  vertical 
arm,  as  well  as  four  1  KW  cartridge  heater  Water  cooling  coils  were 
soldered  onto  each  of  the  sidearms  and  connected  in  series.  The 
sidearms  were  closed  at  the  end  by  pyrex  windows.  The  vertical  arm 
was  connected  to  a  vacuum  pump,  a  pressure  gauge  and  a  supply  of  inert 
gas.  Three-inch  sections  of  wire  mesh  folded  in  cylindrical  shapes 
were  inserted  into  each  of  the  arms  to  serve  as  the  wicks.  Provisions 
were  made  for  several  thermocouples  by  boring  one-inch  deep  holes  in 
the  stainless  steel  cylinder.  The  chromel-alumel  thermocouples  were 
inserted  and  connected  to  calibrated  direct  readout  temperature 
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Figure  9  -  Side  View  of  Crossed  Heatpipe 
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Figure  10  -  Top  View  of  Crossed  Heatpipe 
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controllers.  The  four  1KW  heaters  were  connected  in  a  series-parallel 
arrangement  and  then  wired  to  a  22  amp  Variac.  The  heatpipe  itself 
was  wrapped  in  high  temperature  insulation  and  was  then  enclosed  in  a 
support  capsule,  which  allowed  only  the  arms  to  extend.  Metals  were 
inserted  into  the  heatpipe  from  the  top  by  removing  the  elbow  on  the 
upper  arm  of  the  heatpipe  as  shown  in  Figure  9.  The  heatpipe  was 
operated  at  temperatures  of  600°C  to  1000°C  which  corresponded  to 
Group  I  and  Group  II  vapor  pressures  of  1-200  torr.  The  inert  gas 
pressure  was  varied  from  1-600  torr,  depending  on  the  pressures  of  the 
metals  in  the  heatpipe.  The  inert  gas  pressure  was  always  kept  higher 
than  the  vapor  pressure  of  the  Group  II  metal.  Three  to  four  hours 
were  generally  required  to  bring  the  heatpipe  to  a  stable  temperature 
of  800°C.  Several  days  operation  with  one  fill  of  2-5  grams  of  each 
metal  was  generally  obtained. 

3.  Experimental  Configurations 
a.  Emission  Experiments 

The  purpose  of  the  emissic  i  experiments  was  to  detect  the  presence 
of  LiMg  and  LiCa  excimers  by  observation  of  spontaneous  emission 
resulting  from  the  A2I1  and  B2!  to  X2I  transitions.  To  accomplish 
this  various  laser  sources  were  used  to  excite  the  alkali  metal,  which 
in  turn  was  expected  to  react  with  the  alkaline  earth  producing  the 
excimer . 

The  experimental  configuration  is  shown  in  Figure  11.  A  listing 
of  all  equipment  used  in  all  the  experiments  is  shown  in  Appendix  A. 
The  output  of  an  Ar+  laser,  or  a  dye  laser  pumped  by  an  Ar+  laser  was 
modulated  by  a  mechanical  chopper  at  lOHz-lKHz,  reflected  off  of  two 
turning  mirrors  and  passed  through  a  Brewster  window  into  the  heatpipe 


where  it  was  absorbed  by  the  species  in  the  heatpipe.  The  resulting 
fluorescence  was  focused  by  a  lens,  dispersed  by  a  one-meter  scanning 
monochromator  and  projected  on  a  detector.  A  water-cooled 
photomultiplier  tube  sensitive  from  400-1000nm  and  an  intrinsic  Ge, 
liquid  nitrogen  cooled  solid  state  detector,  sensitive  from  800-1800nm 
were  used  to  observe  the  fluorescence. 

When  the  photomultiplier  was  used,  the  signal  was  first  amplified 
and  then  analyzed  by  a  phase  sensitive  detector.  The  output  of  the  Ge 
detector  was  transmitted  directly  to  the  phase  sensitive  detector. 

The  resulting  output  signals  were  then  recorded  on  the  chart  recorder. 
By  scanning  the  monochromator  at  a  fixed  rate  (A/min)  and  setting  the 
chart  recorder  speed  at  a  fixed  rate  (cm/min),  an  intensity  profile  of 
the  emission  from  the  heatpipe  as  a  function  of  wavelength  was 
obtained. 

Several  considerations  were  important  in  conducting  these  experi¬ 
ments.  First,  the  laser  output  had  to  be  stable  in  both  power  and 
frequency.  Next,  the  f ^  of  the  focusing  lens  system  with  the  heatpipe 
arm  as  its  constraining  aperture  had  to  be  matched  to  the  of  the 
monochromator  for  maximum  light  collection  efficiency.  Additionally, 
the  resolution  of  the  monochromator  was  matched  to  the  product  of  the 
scan  rate  of  the  monochromator  and  the  integration  time  of  the  phase 
sensitive  detector.  Finally,  the  phase  reference  of  the  lockin 
amplifier  was  adjusted  for  maximum  output  signal. 

Several  calibrations  were  performed  to  insure  that  accurate  data 
were  obtained.  The  wavelength  reading  on  the  monochromator  was 
calibrated  using  Ar,  Hg,  and  Ne  lamps  which  emitted  standard  spectral 
lines.  The  monochromator  was  found  to  be  accurate  to  ±1A  over  the 


Figure  11  -  Chemiluminescence  Experiment  Configuration 


range  of  650-1200nm.  The  combined  accuracy  of  the  scan  rate  of  the 
monochrometer (A/min)  and  the  scan  rate  of  the  chart  recorder (cm/min) 
was  determined  to  be  +.01A/cm.  Additionally,  the  combined 
detector/monochrometer  response  curve  was  mapped  from  450nm  to  llOOnm 
using  a  standard  blackbody  source  operating  at  1000°C. 

b.  Absorption  Experiments 

Absorption  experiments  which  utilized  both  a  white  light  source 
and  a  laser  source  were  performed  to  determine  the  presence  of  species 
not  detected  in  emission  and  to  monitor  changes  in  species 
concentrations . 

The  white  light  absorption  experiments  were  performed  by 
modifying  the  setup  in  Figure  11.  The  laser  was  removed  and  a  quartz 
halogen  lamp  was  focused  through  one  arm  of  the  heatpipe.  The 
mechanical  chopper  was  moved  to  a  position  between  the  light  source 
and  the  monochromator.  A  pin-hole  was  also  placed  in  front  of  the 
light  source  to  limit  the  intensity  of  the  source  and  define  the 
source  aperature.  The  same  general  procedures  of  scanning  the  mono¬ 
chromator  as  a  function  of  wavelength  while  recording  the  signal  on  a 
chart  recorder  were  followed.  Observed  variations  in  intensity  as  a 
function  of  wavelength  were  then  correlated  to  absorptions  by  species 
in  the  heatpipe. 

Another  set  of  experiments  were  designed  to  monitor  the 

absorption  of  a  laser  source.  In  these  experiments  an  Ar  laser 

operating  on  a  wavelength  resonant  with  one  of  the  Li^  transitions  as 
46 

shown  in  Table  XI  was  passed  directly  through  the  heatpipe.  Its 


power  was  measured  both  before  and  after  th  heatpipe  at  various 


hie  XI  -  L^tX-B)  Transitions  Resonant  With  Ar  Lines 
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temperatures  with  just  Li  and  with  Li  in  combination  with  Mg.  The 
results  after  measuring  and  accounting  for  losses  on  the  windows 
provided  infomation  on  the  alkali  metal  concentrations  and 
perturbations  to  these  concentrations  due  to  chemical  reactions  or 
large  photon  fluxes. 

D.  Results  and  Analysis 

1 .  Introduction 

Laser  induced  chemiluminescence  experiments  were  initially 
performed  with  many  different  combinations  of  Group  I/II  metals  util¬ 
izing  the  emission  experiment  configuration  previously  described.  Li, 
Na,  and  K  were  tested  in  combination  with  Mg,  Ca ,  Sr,  Ba,  Zn  and  Cd. 

In  all  cases  the  intent  was  to  optically  excite  an  alkali  dimer  or 
interalkali  diatomic  which  would  chemically  react  with  the  Group  II 
metal  to  form  an  excimer  whose  spontaneously  decay  could  be  observed. 
Only  the  combinations  of  Li  and  Mg  and  Li  and  Ca  led  to  the 
observation  of  emissions  that  could  not  be  explained  in  terms  of 
previously  analyzed  alkali  metal  or  alkaline  earth  diatomics.  The 
following  sections  will  describe  the  results  of  the  experiments  with 
Li,  Mg  and  Ca. 

2.  Absorption  Experiments 

An  initial  characterization  of  the  species  present  in  the 
heatpipe  was  performed  to  determine  impurity  levels  prior  to 
conducting  emission  experiments.  Additionally,  interactive  effects 
between  the  Li  and  Mg  or  Ca  such  as  vapor  pressure  suppression  were 
assessed. 

With  the  experimental  apparatus  configured  for  the  white  light 
absorption  experiments,  the  spectrum  shown  in  Figure  12  was  obtained 
with  only  Li  in  the  heatpipe  at  a  temperature  of  900°C  and  432  torr 


total  pressure.  The  published  Li  and  Li2  vapor  pressures  at  900°C  are 

2.5  and  .1  torr,  respectively.  The  absorbing  species  were  identified 
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using  published  spectroscopic  data.  The  primary  absorptions 

resulted  from  the  Li2  (X1Z-A1fl,B1n)  and  Li(2S-2P)  transitions.  There 
were  several  impurities  present  (Ba,  Na,  K) .  Their  concentrations 
were  at  least  two  orders  of  magnitude  less  than  the  Li  concentration 
based  on  comparison  of  the  relative  absorption  areas  for  these  equally 
strong  absorbers.  Trace  amounts  of  Ca  were  also  evident.  Addition  of 
Mg  to  the  heatpipe  resulted  in  the  appearance  of  the  Mg  absorption 
shown  at  457nm.  The  observation  of  Mg  absorption  on  the  forbidden 
lS-3P  transition  was  indicative  of  large  concentrations.  When  either 
Mg  or  Ca  were  added,  no  additional  impurities  were  observed  and  no 
changes  were  noted  in  the  strength  of  the  Li  and  Li2  absorptions. 
Therefore,  there  was  no  apparent  supression  of  the  Li/Li2  concen¬ 
trations  from  the  addition  of  Mg. 

Figure  13  shows  the  results  of  the  laser  absorption  experiments 
for  one  set  of  temperatures.  The  absorption  on  the  Li2(X-B)  transition 
is  linear  in  laser  power  for  the  case  of  Li  alone  and  for  Li  +  Mg 
present  in  the  heatpipe  for  laser  powers  of  .08  watts  to  1.1  watts. 

The  significance  of  the  linear  response  will  be  discussed  in  the  next 
chapter.  Additionally,  there  is  less  than  a  10%  difference  in  the 
laser  absorption  profiles  with  and  without  Mg  present.  Computing  an 
average  value  of  I/Io  for  each  set  of  data  and  applying  Beer's  law, 
[L^Jcrl  varies  by  only  3%.  These  results  indicate  that  the  addition  of 
Mg  does  not  suppress  or  perturb  the  Li^  density  within  the  accuracy  of 
the  experiment  (±10%  due  to  temperature  variations).  The  laser 
absorption  data  supports  the  quantitative  data  from  the  white  light 


44 


INTENSITY  (ARB  UNITS) 


Figure  12(a)  -  Li  Absorption  Spectrum  450  -  520nm 
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Figure  12(b)  -  Li  Absorption  Spectrum  530  -  600nm 
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Figure  12(c)  -  Li  Absorption  Spectrum  610  -  680nm 
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Figure  12(d)  -  Li  Absorption  Spectrum  690  -  760nm 
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Figure  13  *  Laser  Absorption  Data 


absorption  experiments  and  provides  the  additional  information  that  no 
perturbation  is  observed  as  a  result  of  the  absorption  of  the  laser 
light  and  subsequent  excited  state  reactions  involving  Mg.  Therefore 
published  data  on  the  vapor  pressures  can  be  used  to  determine  the 
species  number  densities  at  a  given  temperature. 

2.  Fluorescence  Experiment  Results 
a.  Results 

As  in  the  absorption  experiments,  emission  spectra  were  recorded 

with  only  Li,  Mg,  and  Ca  and  combinations  of  Li  and  Mg  and  Li  and  Ca. 

The  results  of  the  Li  only  experiments  are  shown  in  Figure  14.  The 

spectrum  was  recorded  with  Li  present  at  a  temperature  of  900°C  and 

total  pressure  of  400  torr  with  the  Ar+  laser  operating  on  the  488nm 

line  at  1  watt  output  power.  Various  other  lines  of  the  laser  (476 

and  458nm)  also  resulted  in  observation  of  the  same  emissions  only 

varying  in  intensity.  The  spectra  were  assigned  as  shown  and  the 

observed  emissions  were  identified  as  Li2(B-X),  470-530nm;  Li(2P-2S), 

671nm;  and  Li2(A-X),  630-900nm;  based  on  comparison  of  observed 

43-45 

emissions  to  published  data.  No  emission  was  observed  from  any 

impurities.  Additionally,  no  emission  whatsoever  was  observed  at 
wavelengths  of  930nm-1300nm.  Emission  spectra  taken  under  the  same 
conditions  with  only  Mg  or  Ca  present  showed  no  emission  in  the  range 
of  450nm-1300nm. 

The  combinations  of  Li/Mg  and  Li/Ca  resulted  in  observation 
of  the  Li  emissions  described  above  reduced  in  intensity  50-70%  with 
the  added  features  shown  in  Figures  15  and  16,  respectively. 

Several  experiments  were  performed  to  determine  the  source  of  these 
new  emissions  and  the  mechanisms  by  which  they  were  produced.  In  both 
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Figure  14(a)  -  Li^CB-X)  Fluorescence  460-560nm 
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Figure  15  -  LiMg  Emission 
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Figure  16  -  LiCa  Emission 
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cases  large  amounts  of  impurities  (K,  Na,  Ba)  were  introduced  into  the 
heatpipe.  No  effect  was  observed  on  either  spectrum.  Additionally  Ne 
was  substituted  for  Ar  as  the  inert  gas  again  with  no  effect.  This 
led  to  the  conclusion  that  Li  and  Ca  or  Li  and  Mg  were  the  only 
participants  in  the  emission.  In  addition,  the  dependence  of  the 
intensity  of  the  LiMg  emission  on  laser  power  and  inert  gas  pressure 
was  investigated.  Figures  17  and  18  show  that  the  dependence  on  laser 
power  was  linear  and  there  was  an  inverse  dependence  on  inert  gas 
pressure . 

A  qualitative  investigation  of  the  temperature  dependence  of 
the  emission  intensity  showed  that  the  LiMg  intensity  maximized  at 
900°C  and  was  first  detectable  was  800°C.  The  LiCa  emission  was  first 
visible  at  920°C  and  was  still  increasing  at  the  maximum  temperature 
of  operation  970°C. 

In  addition  to  the  1040nm  spectrum  observed  with  Li  and  Ca 

present  in  the  heatpipe,  the  feature  shown  in  Figure  19  was  also 

observed.  Under  higher  resolution  numerous  bandheads  were  observed. 

The  780nm  emission  exhibited  a  vibrational  spacing  of  ’ =  130/cm, 

u)  "=70/cm  and  AT  =13, 070/cm.  The  observed  emission  did  not  correlate 
e  e 

to  any  predicted  LiCa  transitions.  Additionally,  the  formation  of 

highly  excited  LiCa  would  have  had  to  result  from  energy  pooling 

processes  and  therefore  the  emission  intensity  should  be  quadratic  in 

laser  power.  Figure  20  shows  that  the  power  dependence  is  clearly 

linear.  Vibrational  analysis  suggests  that  the  emission  correlates  to 

a  previously  unobserved  transition  between  the  first  excited  state  of 

Ca9(Ca(3P)+Ca(1S))  since  the  u»  "  agrees  with  previously  determined 

e 
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values  for  the  ground  state  spacing  for  Ca  . 
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17  -  LiMg  Emission  Intensity  Dependence  on  Laser  Power 
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Figure  20  -  LiCa  780nm  Emission  Intensity  Dependence  on  Laser  Power 
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b.  LiMg  and  LiCa  Spectroscopic  Analysis 

Both  the  960nm  LiMg  emission  and  1040nm  LiCa  emission  exhibited 
band  heads  characteristic  of  bound-bound  diatomic  spectra  as  well  as 
continuum  structure  characteristic  of  excimer  emission.  Analysis  of 
the  emission  was  limited  by  several  factors.  Since,  the  excimers  were 
formed  by  chemical  reaction  there  was  no  selective  way  of  producing  a 
given  upper  state  vibrational  level  as  in  a  laser  fluorescence 
experiment  so  that  individual  vibrational  transitions  could  be 
isolated.  Additionally,  the  observed  vibrational  progressions 
included  only  a  few  vibrational  levels.  A  crude  analysis  was  performed 
based  on  the  assumption  that  the  ground  state  is  much  more  weakly 
bound  than  the  excited  state.  Therefore  smaller  spacings  were  assigned 
to  the  lower  state  and  larger  spacings  to  the  excited  states. 
Additionally,  it  was  assumed  that  the  upper  and  lower  state  potential 
curves  were  not  shifted  very  much  from  one  another.  This  results  in 
transitions  such  as  v'=0  to  v"=0,  v'=l  to  v"=l.  If  there  is  a  large 
shift,  the  more  likely  transitions  might  be  v'=4,  5,  6  to  v"=0  for 
example.  The  data  provided  no  way  to  distinguish  between  the  two 
cases.  Using  these  assumptions  the  vibration  assignments  shown  in 
Tables  12  and  13  were  made  for  the  960nm  LiMg  emission  and  1040nm  LiCa 
emission.  The  resulting  spectroscopic  constants  are  shown  in  Tables  14 
and  15  with  the  calculated  constants  from  Tables  1  and  2  in  parenthesis. 

The  experimental  results  for  LiMg  and  LiCa  can  now  be  compared  to 
the  theoretical  predictions  from  Tables  3  through  10  as  shown  in 
Tables  16  and  17.  The  LiMg  results  show  that  the  experimentally  derived 
vibrational  assignment  for  the  upper  state  was  shifted  by  four 


vibrational  levels. 


Table  XII  -  LiMg  Spectral  Assignment (X-nm) 


v' 

v" 

0 

1 

2 


0  1  2 

974.4  951.4 

979.4  955.4 

958.4  935.6 


Table  XIII  -  LiCa  Spectral  Assignments (,\-nm) 
a)  A2n-X2I 


0  1049.9 

1  1056.3 

2  1062.7 

3  1067.6 


b)  B2I-X2I 


0  1047.0 

1  1052.7 

2  1058.7  1038.0 

3  —  1044.1 
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Table  XIV  -  LiMg  Spectroscopic  Constants 


State 

Constant 

X2I 

A2n,B2I 

Tg(ev) 

0 

1.27(1.12,1-56) 

Dg(ev) 

— 

.580  76, .30) 

ui^O/cm) 

50(60) 

250(302,328) 

u,exe^1^cm') 

23(7) 

— 

Table  XV  -  LiCa  Spectroscopic  Constants 


State 

Constant 

T  (ev) 
e 

D  (ev) 
e 

w  (1/cm) 
e 

wexe(1/ cm) 


X2I  A2n  B2I 


0  1.181(1.08)  1.184(1.17) 

.666  .663 

55(93)  —  188097) 

5(7) 
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Table  XVI  -  LiMg  Experimental  vs.  Predicted  Spectrum 


v' 

v" 

Wavelength 

(microns) 

Relative 

Intensity 

A2* 

Xz2 

A  . 
calc 

^obs 

calc 

obs 

3 

0 

1.0060 

— 

— 

1 

1.0011 

---- 

---- 

4 

0 

.9770 

.9744 

.08 

.54 

1 

.9810 

.9794 

.15 

.33 

5 

0 

.9512 

.9514 

.13 

.51 

1 

.9552 

.9559 

.25 

.81 

2 

.9592 

.9584 

.38 

1.00 

6 

0 

.9268 

.38 

— 

1 

.9312 

.54 

— 

2 

.9342 

.9356 

.69 

.22 

7 

0 

.9037 

.77 

— 

1 

.9109 

1.00 

... 
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Table  XVII  -  LiCa  Experimental  vs.  Synthetic  Spectrum 


v' 

v" 

Wavelength 

(microns) 

Relative 

Intensity 

B2I 

X2I 

\  , 
calc 

K  . 
obs 

calc 

obs 

0 

0 

1.0472 

1.0470 

1.00 

1.00 

1 

1 

1.0531 

1.0527 

.53 

.60 

2 

1.0593 

1.0587 

.18 

.  10 

3 

1.0650 

-- 

.03 

-- 

1 

0 

1.0260 

-- 

.  12 

-- 

% 

1 

1.0320 

-- 

.07 

-- 

2 

1.0379 

1.0380 

.34 

.40 

3 

1.0441 

1.0441 

.44 

.50 

A2n 

X2I 

3 

0 

1.0501 

1.0499 

.24 

.30 

1 

1.0569 

1.0563 

.28 

.40 

2 

1.0625 

1.0627 

.24 

.25 

* 

3 

1.0679 

1.0676 

.  18 

.15 

f 

£ 
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Even  after  this  correction,  there  is  poor  agreement  between  the 
intensity  distributions.  The  discrepancies  observed  in  the 
spectroscopic  constants  can  be  clarified  after  matching  the  calculated 
spectra  to  the  experimental  spectra. 

Examining  LiMg  first,  Table  16  shows  the  comparison  of 
experiment  to  theory.  The  first  important  point  is  that  the  observed 
emission  can  be  assigned  to  the  A2n-X2I  transitions.  The  assigned 
transitions,  however,  were  off  by  four  vibrational  levels  in  the  upper 
state  (ie,  the  observed  transitions  are  4,0;  5,0;  etc,  not  0,0;  1,0) 
The  discrepancies  in  the  vibrational  constants  result  from  the 
inability  to  account  for  the  rotational  effects  in  the  experimental 
data.  The  term  in  the  energy  equation  for  a  vibrational  and 
rotational  level  that  includes  coupling  of  vibrational  motion  and 
rotational  motion  is  a  (V+lj)  ( J ( J+l )  )  .  For  large  J  values  (J  *v35 
for  T=900°C)  this  introduces  a  shift  in  the  spectrum  that  is  propor¬ 
tional  to  the  vibrational  quantum  number  which  in  effect  reduces  the 
observed  vibrational  spacing. 

There  are  several  other  apparent  discrepancies  between  the 
theory  and  experiment.  The  theoretical  calculations  indicated  that 
the  most  intense  emission  should  emanate  from  the  LiMg  B2Z  states  in 
the  region  of  700-800nm.  If  the  formation  reaction  is  Li2(B)  and  Mg 
— >  LiMg  and  Li,  then  there  is  11,600/cm  energy  available  for 
population  of  LiMg  states  (Li2  is  bound  by  8900/cm  and  a  488nm  photon 
carries  20,500/cm  energy  resulting  in  11,600/cm  available  energy). 
Since  Tg  of  LiMg  B2Z  is  12,600/cm  this  state  is  not  accessible  to  the 
formation  reaction  and  no  B2I-X2I  emission  should  be  observed. 
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The  reaction  to  form  LiMg(A2Fl)  v'=0  is  2400/cm  exothermic  and 
thermochemically  neutral  for  A2n,  v'=8.  The  exothermicity  of  the 
reaction  may  provide  the  basis  for  explaining  the  apparent  discrepancy 
between  the  observed  and  calculated  A2I1-X2I  intensity  profiles  shown 
in  Table  16.  The  chemical  reaction  may  be  selectively  populating  the 
lower  vibrational  levels. 

An  alternative  explanation  for  the  intensity  variation  is  an 
error  in  the  calculation  of  the  ground  state  potential.  Small  changes 
in  the  shallow  ground  state  potential  will  shift  the  classical  turning 
points  of  the  vibrational  levels  resulting  in  significant  changes  in 
the  Franck-Condon  envelope. 

A  small  shift  of  the  classical  turning  points  of  the  ground 

o 

state  vibrational  levels  to  shorter  internuclear  separation  (.3A) 

o 

either  from  an  inward  shift  of  the  entire  ground  state  well  (3.9A — > 

o 

3-6A;  less  than  a  10%  change)  or  from  a  change  in  shape  would  shift 
the  intensity  maximum  from  the  calculated  v'=9,  v"=0  to  v'=4,  v"=0 
transition. 

The  observed  continuum  emission  from  LiMg  peaks  at  a  wavelength 
of  960-980nm  which  is  consistent  with  the  fact  that  the  most  intense 
bound-bound  emission  results  from  v’=4,  5,  6. 

Excellent  agreement  was  obtained  between  theory  and  experi¬ 
ment  in  the  case  of  LiCa  as  shown  in  Table  17.  Th*  ■•periment  provides 
evidence  that  emission  results  from  population  f  the  A2n  and  B2! 

states.  Since  these  states  lie  lower  in  energy  than  the  corresponding 
states  in  LiMg  (A2n,  Tg  =  8711/cm  and  B2I,  Tg  =  9437/cm)  there  is 
sufficient  energy  available  in  the  formation  reaction  (11,600/cm)  to 
populate  both  states.  Finally,  the  LiCa  ground  state  was  easier  to 


calculate  because  it  has  substantially  more  binding  than  the  LiMg 
ground  state  (580/cm  vs  135/cm)-see  tables  I  and  II.  Both  of  these 
factors  may  have  led  to  closer  agreement  in  the  LiCa  case. 

Figure  21  demonstrates  the  agreement  between  the  calculated  and 
experimental  bound-continuum  emission  spectrum  as  well  as  the 
bound-bound  spectra.  The  continuous  emission  which  originates  from 
transitions  from  the  v'=3,  4,  5  of  the  B2!  state  to  the  continuum, 
peaks  in  the  region  of  lOlOnm  as  predicted  theoretically.  For 
purposes  of  clarity  the  bound-bound  emissions  were  not  weighted  by  the 
large  number  of  rotational  transitions  emitting  times  their  Doppler 
width.  This  results  in  the  apparent  amplitude  discrepancy  between  the 
bound-bound  and  bound- continuum  theoretical  spectra. 

Finally,  it  should  be  noted  that  there  was  good  agreement  between 
the  observed  and  calculated  spectroscopic  constants  with  the  exception 
of  the  ground  state  As  in  the  case  of  LiMg,  this  discrepancy  can 

be  accounted  for  by  consideration  of  the  large  vibrational-rotational 
coupling  constant  for  the  ground  state,  together  with  the  larger 
rotational  quantum  numbers  (Jmax=^0)  resulting  from  the  high 
temperature. 

E.  Summary  of  Results 

A  study  of  combinations  of  Group  IA/IIA  metals  resulted  in  the 
identification  of  two  new  emitters,  LiMg  and  LiCa.  These  species  were 
produced  as  a  result  of  an  initial  excitation  of  Li2(B)  and  subsequent 
chemical  reaction.  The  observed  emissions  were  correlated  with  the 
results  of  ab  initio  calculations  with  good  agreement.  Several  other 
pieces  of  data  were  obtained  on  LiMg  which  will  affect  the  kinetics 
analysis  as  explained  in  the  next  chapter. 
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Figure  21  -  LiCa  Theoretical  Spectrum  Compared  to  Experimental  Results 
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Chapter  III  -  LiMg  Kinetics 


A.  Introduction 

This  chapter  presents  an  analysis  of  the  kinetic  processes 

relevant  to  gas  phase  Li/Li^  and  Li/L^/Mg  systems.  A  kinetic  model 

is  developed  based  on  the  experimental  observations  from  the  previous 

chapter  and  published  studies  of  the  equivalent  Na/Na^  and  K/K^ 

48 

systems.  Time  dependent  and  steady  state  analyses  of  the  rate 
equations  are  presented  and  correlated  with  experimental  measurements 
of  time  delays  and  relative  intensities  to  determine  the  rate 
constants . 

B.  Theory 

1.  Species  Concentrations  in  a  Li/Mg  gas. 

The  primary  constituents  of  an  equilibrium  gas  phase  mixture 
of  lithium  and  magnesium  are  Li,  Mg,  Li2,  Mg2  and  LiMg.  The  contributions 
of  polyatomics  (Li^,  Li^Mg,  etc.)  to  the  total  pressure  are  assumed  to 
be  negligible^ since  not  only  is  their  binding  energy  expected  to 
be  small  but  their  formation  rates  are  also  expected  to  be  small  since 
the  potential  reactants  (LiMg,  Li^,  etc.)  are  expected  to  be  present 
in  small  concentrations.  The  gas  phase  concentrations  of  Li  atoms,  Li 
dimers  and  Mg  atoms  have  been  measured  to  within  15%  uncertainty  in 
pure  Li  or  pure  Mg  gases  from  200°C  to  1000°C  as  described  in  Nesmeyanov~> * 
and  are  best  characterized  by  the  expressions  shown  in  Appendix  B. 

Based  on  the  experimental  data  shown  in  Chapter  2,  it  will  be  assumed 
throughout  this  study  that  the  Li2,  Li  and  Mg  densities  can  be  determined 
directly  from  Nesmeyanov's  expressions  for  their  unperturbed  equilibrium 
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vapor  pressure  curves. 


There  is  no  published  data  on  the  equilibrium  concentrations  of 

LiMg  and  minimal  data  on  Mg2 .  Mg2  has  been  observed  in  a  pure  Mg  gas 

52 

and  its  ground  state  was  determined  to  be  bound  by  400/cm. 

Utilizing  this  information,  the  Mg2  density  in  a  pure  Mg  gas  was 

calculated  as  shown  in  Appendix  C.  The  density  ratio  Mg2/Mg  was 

-4 

determined  to  be  1.6  X  10  at  900°C  which  was  a  typical  experimental 
operating  temperature.  Equilibrium  arguments  can  be  used  to  estimate 
the  effect  of  Li  on  Mg2  densities.  The  LiMg  equilibrium  concentration 
can  be  obtained  similarly.  The  Mg2  density  is  coupled  to  the  Li2 
density  through  LiMg  as  shown  below: 

Mg2  +  Li  ■*-»  LiMg  +  Mg 
Li^  +  Mg  LiMg  +  Li 

The  first  reaction  is  nearly  thermoneutral  since  the  equilibrium 

constant  is  exp(-AE/kT)  =  .3  for  kT  =  800/cm  at  900°C  and  AE  is 

determined  from  the  difference  of  the  Mg2(400/cm)  and  LiMg(135/cm) 

bond  energies.  Therefore  nearly  equal  concentrations  of  Mg^  and  LiMg 

will  be  maintained.  The  second  reaction  is  strongly  driven  to  the 

left  since  Li2  is  bound  by  8900/cm  and  the  resulting  equilibrium 

constant  is  2.1  x  10  \  In  effect  the  Li2  concentration  suppresses 

the  LiMg  concentration  which  in  turn  keeps  the  Mg2  concentration  low. 

This  effect  was  quantified  by  application  of  the  Gibbs-Duhem 
53 

relationship  as  shown  in  Appendix  D.  The  results  were  Mg2/Mg  and 
LiMg/Mg  ratios  of  1.2  x  10  ^  and  2  x  10  ^  at  the  reference  temperature 
of  900°C.  Therefore,  the  result  of  adding  Li  to  a  pure  Mg  gas  is  to 
suppress  the  Mg2  density  by  an  order  of  magnitude  and  also  produce  a 
small  equilibrium  concentration  of  LiMg.  Table  XVIII  summarizes  the 
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absolute  magnitude  of  all  the  species  concentrations  in  the  mixed  LiMg 
gas.  The  knowledge  of  these  relative  species  concentrations  will 
provide  the  basis  for  determining  the  importance  of  competing 
processes  in  the  optically  excited  Li/Mg  system. 


Table  XVIII  -  LiMg  Gas  Phase  Concentrations  at  900°C 


Species 

Concentration 

Li 

9.0 

X 

10^/cm' 

Li2 

4.8 

X 

1015/ cm 

LiMg 

1.8 

X 

10*  Vcm 

Mg 

9.0 

X 

10^/cm 

Mg2 

1.1 

X 

10^/cm 

2.  Li2(X)  Perturbation  Due  to  Photon  Flux. 

The  data  presented  in  Chapter  2,  Figure  14,  demonstrated  that  the 
LigCX)  density  at  900°C  was  unperturbed  (±1 0%)  by  optical  excitation 
from  a  488nm,  1  watt  Ar+  laser  beam.  A  theoretical  estimate  of  the 
magnitu  '  ■  of  the  perturbation  can  be  obtained  from  a  rate  equation 
analysis  uf  the  Li2(B)  concentration  including  both  optical  excitation 
and  decay  terms  while  ignoring  stimulated  emission: 

.l|Li  fB)]/dt  =  Io[  Li  ^  (X )  )/hv  +  lLi,(B)]/T 


2.1 


where  I  is  the  intensity  of  the  laser  (]W/cm2),  O  is  the  effective 
absorpti?  cross  section  (cm2)  of  the  Li^CX-B)  transition,  hv  is  the 
photon  energy  (20500/cm),  and  1  / X  is  the  L^tB)  decay  rate  due  to  both 
radiative  and  collisional  processes.  Neglecting  the  collisional 

g 

effects,  l/l  is  shown  in  Appendix  E  to  be  10  /s.  The  product  alLi^fX)] 
was  experimentally  determined  in  Chapter  2  to  be  .12/cm  at  900 °C. 

Using  these  results  in  equation  1  yields  [ Li^ (B) )  =  3  xlO  /cm  and 
(Li^fB)  ]/ [^^(X) )  =  10  ^  where  the  density  of  L^tX)  corresponds  to 
the  unperturbed  equilibrium  vapor  pressure.  Other  terms  in  the  rate 
equation  which  account  for  reactive  losses  of  Li2(B)  such  as 
collisions  with  Li  atoms,  will  further  suppress  the  excited  dimer 
density.  The  (L^^)]  density  estimate  derived  here  is  therefore  an 
upper  limit  which  shows  that  the  perturbation  of  the  ground  state  due 
to  the  optical  excitation  does  not  exceed  one  part  in  10^. 

3.  Reactions  Involving  Excited  State  Species 

The  excited  Li2(B)  molecules  produced  by  laser  excitation  may 

transfer  energy  to  or  chemically  react  with  other  species  as  a  result 

of  collisions.  Such  processes  have  been  observed  in  both  the  Na/Na,, 

2 

and  K/I<2  systems.  In  the  Na/Na^  system,  Na2(B)  reacts  with  Na(  S)  to 
produce  Na2(X)  and  Na(2P).^’^  Additionally,  the  reaction  of  Na2(X) 
with  Na(2P)  has  been  shown  to  form  Na(2S)  and  Na2  (a,A).^  By 
analogy,  these  same  processes  are  postulated  to  occur  in  the  Li/Li2 
system  since  the  energy  level  structures  and  symmetries  of  the  systems 
are  the  same  with  one  exception.  As  shown  in  Figure  22  the  dissociative 
3I  state  in  predissociates  the  a3Il  state  at  very  low  vibrational 
levels  which  does  not  happen  in  Na^.  Since  predissociation  rates  are 
much  faster  than  collisional  or  radiative  rates,  processes  involving 


72 


3 

Li^(a  n )  as  a  reactant  can  be  ignored  in  the  Li/Li^  system.  The 
important  reactions  for  the  Li  system  are  shown  below. 


Excitation:  Li2(X)  +  hv  -*  Li2(B)  (1) 

Reaction:  Li2(B)  +  Li(2S)  ->•  Li2(X)  +  Li(2P)  (2) 

Li (2P)  +  Li2(X)  ->  Li2(A,a)  +  Li  (2S)  (3) 

Li  (2P)  +  Li  (2S)  -*  Li  (2S)  +  Li(2S)  (4) 

Radiation:  Li2(B)  •*  Li2(X)  +  hv  (5) 

Li2(B)  -*•  Li2(A)  +  hv  (6) 

Li2(A)  -»  Li2  (X)  +  hv  (7) 

Li  (2P)  -»  Li  (2S)  +  hv  (8) 


These  reactions  account  for  all  the  emissions  that  were  observed  in 
the  fluorescence  of  Li  vapor  and  reported  in  Chapter  2. 

Since  the  experiments  in  this  thesis  were  performed  in  an 
environment  which  included  an  inert  gas,  the  inert  gas  deactivation 
reactions  must  also  be  considered  in  the  specification  of  the  kinetic 
model  as  shown  below: 

Li  (B)  +  Ar  -»  Li2(X)  +  Ar  (9) 

Li  (A)  +  Ar  ->  Li2(X)  +  Ar  (10) 

Li  (2P)  +  Ar  ->  Li  (2S)  +  Ar  (11) 

Electronic  to  translational  energy  conversion  rates  are  generally 
small  (10  ^  to  10  ^cm^/s)  and  therefore  these  reactions  should  have  a 
small  effect  on  the  overall  kinetics. 
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The  model  that  has  been  developed  for  the  Li  system  ignores 

reactions  involving  collisions  of  two  excited  states  species  to 

form  a  doubly  excited  species  (energy  pooling) and  three  body 

reactions,  as  does  the  previous  work  on  K  and  Na  which  was  conducted 

in  the  same  density  regimes  as  the  experimental  work  in  this  thesis. 

There  are  several  Li/Li^  states  accessible  to  energy  pooling  reactions 

2 

of  Li^CB),  L^(A)  and  Li (  P) .  None  of  these  were  observed  in 

emission.  This  result  is  consistent  with  expectations  in  that  gas 

-10  3 

kinetic  formation  rates  (10  cm  /s)  limit  the  pooling  reactions 

9  3  10 

involving  excited  states  species  with  densities  of  10  /cm  -  10 

3 

/cm  to  form  doubly  excited  states  at  rates  of  .1/s  -  1/s  while  the 
rates  for  forming  singly  excited  species  are  typically  10~*  -  10^/s. 

Three  body  recombination  reactions  such  as  Li  +  Li*  +  M 
Li^  (A)  +  M  are  also  ignored.  At  typical  total  pressures  of  100  torr 
the  three  body  rates  are  several  orders  of  magnitude  less  than  the 

corresponding  two  body  rates,  since  the  product  of  a  typical  density 

1 S  3  32  6 

(10  /cm  )  and  a  typical  three  body  rate  of  10  cm  /s  is  much  less 

-12  3  -10  3 

than  the  typical  two  body  rate  of  10  cm  /s  to  10  cm  / s . 

The  addition  of  Mg  to  an  optically  excited  Li  system,  creates 
several  other  possible  kinetic  channels  including  those  which  form  the 
excimer  LiMg(A).  The  possible  new  reactions  are  listed  below: 


Reaction:  Li2(B)  +  Mg(*S)  ■*  Li2(X)  +  Mg  (3P)  (12) 

L i 2 ( B )  +  Mg(JS)  -  LiMg(A)  +  Li(zS)  (13) 

Mg2(1I)  +  Li(2P)  -►  LiMg(A)  +  Mg(2S)  (14) 

Mg  (3P)  +  Liz(X)  -  LiMg(A)  +  Mg(1S)  (15) 

LiMg(A)  +  Li (ZS)  -  Li2(A,a)  +  Mg(1S)  (16) 

Radiation:  LiMg(A)  -*  LiMg(X)  +  hv  (17) 
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Quenching:  Li2(B)  +  Mg(xS)  ->  Li2(A,a,X)  +  Mg^S)  (18) 

LiMg(A)  +  Ar  -»  LiMg(X)  +  Ar  (19) 

Dissociation:  LiMg(A)  +  Q  -*  Li  +  Mg  (20) 

These  reactions  follow  from  consideration  of  reaction  exoergicity, 
as  well  as  the  reaction  correlation  diagram  shown  in  Figure  23.  Several 
of  the  above  reactions  can  be  shown  to  be  negligibly  important  in  the 
specification  of  the  LiMg  kinetic  model.  Reaction  (12),  the  production 
of  Mg(3P)  from  collisions  of  Mg  with  Li2(B),  is  3000/cm  endothermic 
(kT  is  815/cm  at  900°C)  and  spin  forbidden.  Additionally,  this 
reaction  competes  with  the  reaction  forming  the  excimer  (13)  which  is 
exothermic  and  spin  allowed.  Also,  there  was  no  observation  of 
Mg(3P)  emission  in  the  experiments.  For  these  reasons,  reactions 
involving  Mg(3P)  (12,15)  will  be  ignored.  Reaction  (14)  involving 
formation  of  LiMg(A)  by  collisions  of  Li(2P)  with  Mgi?  (X 1  i )  can  be 
similarly  ignored,  since  the  reaction  is  spin  forbidden  and  the  small 
densities  of  both  Li(2P)  and  Mg^fX1!)  present  will  combine  to  make  this 
reaction  many  orders  of  magnitude  smaller  than  the  reactions  forming 
LiMg(A)  by  collisions  of  Mg  with  Li2(B)  as  shown  below. 

A  steady  state  analysis  of  the  production  of  Li^P)  through  col¬ 
lisions  of  Li(2S)  with  Li2(B)  yields  the  following  expression  lor 
Li (2P) : 

[Li (2P) ]  =  [ Li 2 ( B ) ) [Li ]Kf/A(Li ) 

where  is  the  formation  rate  and  A(Li)is  the  radiation  trapped  decay 
rate.  The  maximum  value  of  [Li2(B)]  was  previously  determined  to  be  3  x 
109/cm3,  (Li]  is  9  x  1016/cm3  at  900°C  from  Table  XVIII  and  A(Li)  is  1 
x  10^/s  by  analogy  with  Na.  By  assuming  a  gas  kinetic  formation  rate  of 
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I MTERNUC LEAR  DISTANCE 


Figure  23  -  LiMg  Reaction  Symmetries 
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10  ^cm3/s,  it  is  found  that  an  upper  limit  to  [ Li ( 2 P ) ]  is  3  x 
lO^/cm3  .  With  this  information  the  relative  formation  rates  of 
LiMg(A)  by  reactions  (13)  and  (14)  can  then  be  compared.  Therefore 
[Li (2P) ] [Mg2 ] K i / [ Li2 (B) ] [Mg]K2  equals  10  ^  assuming  Kj  =  K^. 
Additionally  reaction  (14)  should  also  produce  LiMg(B)  emission  which 
is  not  observed.  Therefore  reaction  (14)  will  be  ignored  in  the 
remainder  of  the  analysis. 

Reaction  (18),  the  direct  quenching  of  Li2(B)  to  Li2(X)  is 
unlikely  since  there  are  no  intermolecular  curve  crossings  connecting 
the  states.  Quenching  of  L^tB)  to  the  a3n  state  is  spin  forbidden. 
Quenching  to  the  A 1 n  state  may  be  accomplished  through  curve  crossings 
in  the  triatomic  formed  in  the  collision  of  Li2(B)  with  Mg.  Since  the 
LiMg  excited  states  lie  intermediate  to  the  Li2(B)  and  Li2(A)  states 
there  is  also  a  direct  route  along  intermolecular  curves  by  which 
Li 2 (B )  can  be  converted  in  Li2(A)  as  shown  in  Figure  16.  This 
reaction  is  therefore  likely  to  have  a  significant  rate  in  the  LiMg 
kinetic  mode  1 . 

Finally,  the  possibility  of  collisional  dissociation  of  LiMg(A)  - 

reaction(20)  must  be  considered.  Since  LiMg(A)  is  bound  by  only 

4000/cm  and  has  a  lifetime  of  several  p  seconds  (Appendix  E),  at  a 

typical  pressure  of  100  torr  it  may  undergo  numerous  collisions 
6 

(10  col 1/torr-sec  x  102  torr  x  lpsec  =  100  collisions)  prior  to 
radiative  decay.  Since  5  consecutive  energy  addition  collisions  would 
dissociate  LiMg(A)  for  kT  of  800/cm  reaction  (20)  must  be 
considered.  The  same  reactions  are  ignored  for  Li2(B)  and  Li2(A) 
because  their  radiative  lifetimes  are  two  orders  of  magnitude  shorter 
than  LiMg(A)  which  allows  at  most  only  one  collision  on  the  average 
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before  radiation.  If  the  lifetime  of  the  LiMg(A)  state  is  not 
determined  by  radiative  decay  but  by  reactive  or  collisional  decay 
which  significantly  shortens  the  lifetime  then  the  thermal 
dissociation  process  may  be  ignored  in  LiMg  also. 

As  in  the  case  of  the  mixed  Li  and  Ar  system,  the  effects  of 
three  body  reactions  were  considered  negligible.  The  experimental 
data  shown  in  Figure  18  demonstrated  that  the  LiMg(A)  emission 
intensity  increased  with  decreasing  pressure  in  the  range  of  50-425 
torr.  If  a  three  body  reaction  was  responsible  for  the  excimer 
formation,  there  would  have  been  at  least  an  initial  increase  in 
intensity  with  increasing  pressure  followed  by  a  decrease  when  the 
inert  gas  quenching  rate  equaled  the  product  of  the  third  body  density 
and  the  three  body  formation  rate. 

The  possibility  of  energy  pooling  reactions  in  the  Li/Mg  system 
is  also  ignored  for  the  same  reasons  as  in  the  pure  Li  system;  the 
production  rates  of  doubly  excited  species  are  extremely  small  due  to 
small  (108/cm3  -  10lo/cm3)  excited  state  densities.  Additionally 
there  was  no  observation  of  emission  from  doubly  excited  species. 
Initially,  the  possibility  existed  that  the  LiMg  emission  occurred  as 
a  result  of  an  energy  pooling  process.  However,  the  linear  power 
dependence  of  the  emission  intensity  shown  in  Figure  17  contradicts 
this  hypothesis  since  emission  resulting  from  a  species  produced  as  a 
result  of  energy  pooling  would  have  shown  a  quadratic  dependence  on 
the  laser  power. 

In  addition  to  considering  the  reactions  involved  in  the  Li/Mg 
system,  the  effects  of  radiation  trapping  and  diffusion  need  to  be 
examined.  In  the  thesis  experiments,  a  cyclindrical  volume  of  Li/Mg 
gas  was  excited  along  the  axis  by  a  laser  beam.  The  emission  from 
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this  cylinder  was  viewed  from  a  direction  perpendicular  to  the  axis. 

The  emitted  radiation  therefore  had  to  pass  through  several  centi- 

17  3 

meters  of  the  Li/Mg  gas.  The  typical  gas  densities  at  900°C  were  10  /cm 
163  183 

Li,  10  /cm  Lij  and  10  /cm  Mg.  Since  radiation  trapping  is  known  to 

13  3 

have  an  important  effect  on  atomic  emitters  at  densities  of  10  /cm 
the  Li(2P-2S)  radiation  was  expected  to  be  significantly  radiation 
trapped.  For  the  same  transition  in  Na  the  trapped  lifetime  has  been 
calculated  to  increase  from  a  natural  lifetime  of  16nsec  to  a  constant 
value  of  lOpsec  over  the  pressure  range  of  . 1  -  10  torr."^  In  this 
same  temperature  range,  Na2  radiation  is  not  significantly  radiation 
trapped  since  the  dimer  density  is  20-100  times  lower  than  the  atomic 
density  and  the  number  of  potential  absorbers  per  unit  frequency  is 
substantially  decreased  due  to  rotational  and  vibrational  dilutions. 

These  same  assumptions  will  be  applied  to  the  Li/Mg  system.  The 
lifetimes  of  molecular  emitters  will  be  assumed  unperturbed  by 
radiation  trapping  while  the  trapped  lifetime  of  Li(zP)  will  be 
measured  in  the  experiments. 

It  has  been  shown  that  optical  pumping  in  the  Li/Li^  system  does 
not  significantly  perturb  the  concentration  of  Li^-  However,  when  Mg 
is  present  an  additional  complication  occurs.  Whereas  in  a  pure  Li 
system  the  optically  excited  Li^CB)  molecules  eventually  decay  back  to 
L^CX),  with  Mg  present  the  dimers  may  be  dissociated  due  to  the 
excimer  formation  reactions.  Therefore  in  steady  state  in  the  optically 
pumped  region  the  dimer  concentration  may  be  suppressed  to  the  point 
where  three  body  recombination  and  diffusion  balance  the  loss  of 
dimers  caused  by  excimer  formation.  Since  diffusion  and  recombination 
are  normally  slow  processes  there  is  a  possibility  that  the  Li^ 
concentration  will  be  substantially  depressed.  Figure  13  demonstrated 
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that  experimentally  this  was  not  observed.  It  is  difficult  to  analyze 

this  problem  a  priori  since  the  rate  of  dimer  loss  due  to  excimer 

formation  is  not  known  and  it  is  difficult  to  assess  the  effects  of 

reactions  (1)  and  (2)  which  tend  to  enhance  the  rate  of  formation  of 

Li^  in  addition  to  suppressing  Mg^ . 

Radial  diffusion  of  ground  state  species  is  negligible  since  the 

ground  state  density  gradients  are  small  as  was  previously  shown.  The 

excited  species  density  gradients  however  are  much  larger  and  since 

they  are  produced  in  the  center  region  they  will  diffuse  out. 

”6  *  8 

However,  the  short  lifetimes  of  these  species  (10  s  -  10  s)  will 

prevent  them  from  diffusing  very  far  before  radiative  decay.  Assuming 

2 

a  hard  sphere  model  the  diffusion  rate  D  is  just  v/2[N]nr  where  v  is 
the  thermal  velocity  v  =  -JlkT/m,  { N]  is  the  number  density  and  r  the 

o 

species  radius.  Using  typical  values  for  the  100  torr,  900  C  regime  D 

2  *6 
=  1  cm  /s.  An  excited  state  species  with  a  lifetime  of  T  <  10  sec 

diffuses  typically  -JY)X=  10  cm  prior  to  decay.  Therefore  since  the 

laser  beam  is  typically  10  ^cm  in  diameter  diffusive  loss  of  the 

excited  species  can  be  ignored. 

The  following  figures  summarize  the  important  features  of  the 

LiMg  kinetic  model.  Figure  24  shows  the  energy  flow  in  an  optically 

excited  Li/Mg  system.  Figure  25  summarizes  the  important  reactions. 

The  approximations  that  have  been  made  in  defining  the  model  are 

listed  in  Figure  26.  The  rate  equations  corresponding  to  the 

reactions  listed  in  Figure  25  are  summarized  in  Figure  27. 

4.  Steady  State  Analysis 

Using  the  rate  equations  from  Figure  20  steady  state  solutions 
were  determined  for  the  various  excited  state  emitters  for  the  cases 
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of  a  pure  Li  system  and  a  mixed  Li/Mg  system  as  shown  in  Figure  21.  In 
these  so’  ions  the  ground  state  species  densities  are  assumed  to  be 
the  equilibrium  concentrations  specified  by  the  temperature.  These 
solutions  will  be  related  to  the  results  of  experiments  where  the 
relative  emission  intensities  of  L^CA),  Li^CB)  were  measured  in  pure 
Li  and  mixed  Li/Mg  systems.  The  measured  intensity  ratios  were  used 
to  infer  the  equivalent  concentration  ratios.  Using  this  experimental 
data  several  of  the  rate  constants  in  the  equations  in  Figure  21  were 
determined . 
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Reaction 


Rate  Constant  Designation 


Excitation 

Li2(X)+hu^Li2(B) 

Primary 

Li2(B)+Mg^LiMg(A)+Li 

Kj (cm3/s) 

Li2(B)+Li^Li2(X)+Li(2P) 

K2(cm3/s) 

Secondary 

Li(2P)+Li2(X)^Li2(A,a)+Li 

K3 ( cm3 / s ) 

LiMg(A)+Li-»Li2(A,a)+Mg 

K4(cm3/s) 

LiMg(A)+Ar-»Li  (2P)+Mg+Ar 

K5  (  CfTl3/ S  ) 

Radiation 

Li  (2P)-»Li+hv 

Aj (1/s) 

Li2(B)^Li2(X)+hv 

A2(1/s) 

Li2(A)-^Li2(X)+hv 

A3 ( 1 / s ) 

Li2(B)->Li2(A)+hV 

A4(1/s) 

LiMg(A)-»LiMg+hv 

A5(1/s) 

Quenching 

Li2 (B)+Ar^Li2 (X)+Ar 

K6 

Li2(A)+Ar->Li2(X)+Ar 

K7 

LiMg(A)+Ar->LiMg(X)+Ar 

K8 

Li  (2P)+Ar  ,Li->Li+Ar  ,  Li 

K9 

Li2(B)+Mg^Li2(A)+Mg 

K10 

*Na  System  Rates 


Figure  25  -  LiMg  Reactions 
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Literature  Value 

None 

(2xl0”9 )* 

(3 . 4xl0~9 ) '•' 

None 

None 

3 . 72xl07 

lxlO^(Appendix  E) 
5.56xl07 

lxlO7 (Appendix  E) 

1 . lxlO^(Appendix  G) 

None 

None 

None 

None 

None 


1.  Ground  State  Densities  Li,  Li^,  Mg  are  Unperturbed  by: 

a)  Li-Mg  Interactions 

b)  Laser  Interactions 

2.  Pooling  Reactions  are  Insignificant 

3.  Three  Body  Reactions  are  Insignificant 

4.  Li/Li^  Kinetics  are  Analogous  to  Na/Na^  Kinetics 

5.  Endothermic  Reactions  Considered  Insignificant 

6.  Diffusive  Effects  are  Unimportant 

7.  Radiation  Trapping  Affects  Only  Atomic  Emitters 

8.  Mg2>  Li^Ca)  Reactions  Insignificant 

Figure  26  -  List  of  Assumptions 
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d[Li2(A)]/dt  =  [LirP)][Li2(X)]K3  +  (LiMg(A) ]  (Li]K^ 

+  [Li2(B))AA  -  [Li2(A)](A3  +  K? [Ar] ) 

d[Li2(B)]/dt  =  Io[Li2(X)]/hv  -  [Li2(B)](A2  +  A4  +  [Mg] 
+  K2[Li]  +  K6  [Ar)) 

d [Li (2P) ] /dt  =  [Li2(B)][Li]K2  -  [Li (2P) ]( [Li2(X) ]K3  +  A] 
+  [Ar)Kg)  +  [LiMg(A)]K5[Ar] 

d[LiMg(A)]/dt  =  [Li2(B))[Mg]K1  -  [ LiMg(A) ] ( [Li ]K^  +  A, 

+  [Ar] (Kg  +  K5) 


Figure  27  -  LiMg  System  Excited  State  Rate  Equations 


5.  Time  Dependent  Analysis  of  Li/Mg  Kinetic  System 

If  _  ■  laser  excitation  source  which  produces  L^CB)  is  modulated 

at  a  frequency  iu  then  the  resulting  products  of  the  excitation  will 

also  exhibit  a  similar  periodic  time  dependent  behavior.  Since  the 

reaction  products  in  the  system  are  formed  through  collisional 

processes  and  decay  by  both  radiative  and  collisional  processes  taking 

finite  amounts  of  time,  the  system  response  will  lag  the  excitation. 

These  time  lags  or  phase  shifts  can  be  measured  and  related  to  kinetic 

and  radiative  rates.  The  technique  of  time  delay  or  phase  shift 

measurement  has  been  used  extensively  for  the  determination  of  the 

55-58 

lifetimes  of  excited  states.  In  many  cases  of  experimental 

interest  a  simple  relationship  exists  for  the  phase  shift  (|)  of  the 
emitter  relative  to  the  exciter:  tan<[>  =  U)T  where  w  is  the  modulation 
frequency  and  T  the  excited  state  lifetime.  In  general,  however,  the 
response  of  a  kinetic  system  involving  many  kinetically  active  species 
will  not  be  as  simple.  In  fact,  it  may  be  too  complex  to  be  of  use  in 
determining  kinetic  rates  except  in  certain  limiting  cases  wich 
simplify  the  input/output  relationship.  The  simplest  case  occurs  when 
the  excitation  or  perturbations  are  sufficiently  small  that  the  system 
responds  in  a  linear  fashion.  In  order  to  justify  the  use  of  linear 
response  or  small  signal  theory  to  characterize  the  response  of  the 
LiMg  system,  it  is  necessary  to  establish  that  the  magnitude  of  the 
nonlinear  terms  in  the  rate  equation  are  small  when  compared  to  the 
magnitude  of  the  corresponding  linear  terms. 

The  previous  analysis  and  experimental  measurements  of  the  Li/Mg 
system  have  demonstrated  that  many  of  the  nonlinear  processes  are 
unimportant  in  modeling  the  LiMg  system.  Experimentally  the  absorption 
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of  laser  power  was  observed  to  be  linear  in  laser  intensity  both  with 
and  without  Mg  present.  Additionally  the  LiMg  emission  intensity 
showed  a  linear  dependence  on  laser  power.  Further,  the  analysis  in 
the  previous  several  sections  supported  these  experimental 
observations  and  also  showed  that  the  nonlinear  processes  of 
diffusion,  energy  pooling,  and  three  body  reactions  can  be  ignored. 
However,  even  with  the  elimination  of  these  processes,  there  remain 
certain  nonlinear  terms  in  the  rate  equations  shown  in  Figure  27.  The 
nonlinearities  arise  from  terms  which  contain  products  of  excited 
state  and  ground  state  densities  as  well  as  products  of  ground  state 
densities  and  the  laser  intensity.  The  ground  state(G),  excited 
state(E),  and  laser  intensity(I)  terms  are  all  functions  of  time  that 
have  a  constant  term  and  a  time  varying  part  as  shown  below.  The 
magnitude  of  the  time  varying  portions  are  always  less  than  the 
constant  terms  since  G,E  and  I  are  positive  definite  quantities. 

G(t)  =  Gq  +  AG(t) 

E(t)  =  EQ  +  AE(t) 

I (t)  =  IQ  +  AI(t) 

The  nonlinear  terms  in  the  rate  equations  which  arise  from  products  of 
these  expressions  are  AG(t)AE(t)  and  AG(t)AI(t).  Analysis  in  the 
previous  sections  has  shown  that  typically  AG/Gq  was  approximately  10 
and  AG  is  approximately  equal  AE .  Since  the  laser  was  at  least  70% 
modulated  AI  ~  1^.  Therefore,  all  terms  such  as  AG(t)AE(t)  and 
AG(t)AI(t)  are  always  small  compared  to  terms  like  AE(OGq  and 
AI(t)G0<  The  constant  terms  G^Eq,  etc.,  were  excluded  experimentally 
by  accepting  only  the  time  varying  portion  of  the  output.  The  cross 
terms  AG(t)EQ  and  AG(t)lQ  can  also  be  excluded  when  compared  to  terms 
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like  AE(t)GQ  or  AI(t)GQ  since  AE(t)/EQ  and  Al(t)/IQ  =1  and  G^  >>  Eq. 

Applying  these  concepts  to  the  rate  equations  in  Figure  20 
simplifies  the  nonlinear  differential  equations  to  a  series  of  coupled 
linear  differential  equations.  If  each  species  density  and  the  laser 
intensity  is  then  expanded  in  a  Fourier  series  and  is  substituted  into 
the  rate  equations,  the  differential  equations  are  reduced  to  complex 
algebraic  expressions.  This  method  of  solution  is  applicable  because 
of  the  superposition  principle  which  applies  to  linear  systems  and  the 
orthogonality  of  the  basis  functions  which  make  up  the  Fourier 
expansion.  Then  examining  only  those  terms  at  the  fundamental 
frequency,  applying  the  previous  analysis  regarding  elimination  of 
negligibly  small  terms  and  using  the  Euler  identity,  an  expression  for 
the  phase  shift  (j)  can  be  obtained.  In  the  case  where  the  system  is 
first  order  then  the  simple  result  tan  <])  =  lot  is  obtained  where  T 
represents  the  lifetime  of  an  excited  state.  These  procedures  can 
further  be  clarified  by  using  the  L^OO  rate  equation  as  an  example. 

The  rate  equation  for  Li^CB)  can  be  solved  by  Fourier 
transformation  of  the  periodically  modulated  laser  and  all  the  species 
densities  as  shown  below: 

I  (t)  =  IQ  +  Iiel(uit)  +  .  .  . 

[Li2(B)j  =  [Li2(B)]Q  +  [Li2(B)]1  ei(U,t+4>B)  +  *  *  * 

[Li]  =  [ Li ] o  +  [Li]  j  ei(u,t+*L)  +  '  *  • 

[Li2(X)l  =  [Li2(X)]Q  +  [Li2(X)]1  ‘ 

The  phase  of  the  laser  is  arbitrarily  set  to  zero  while  all  other 


phases  are  defined  to  be  phases  relative  to  the  laser.  All  the 


coefficients  are  real  time  independent  quantities.  By  excluding  the 
nonlinear  terms,  cancelling  e1UJt  in  each  term,  and  keeping  only  the 
terms  at  frequency  ut,  as  previously  discussed,  the  rate  equation  for 
the  Li2(B)  density  converts  to  an  algebraic  expression: 

iw  [Li2(B)]1ei(,,)B)  =  Io(o/hv)  [ Li2 (X) )  2 e 1  (4>X } 

+  I1(o/hv)[Li2(X)]oei(<{,La)  +  [Li2(B)J, 

(A2  +  A4  +  K2  [ Li ) o  +  K6  [Ar]o)el(<|,B) 

+  lLi2(B)}Q  [Li]1K2ei(<^L) 

The  term  I q [ Li2 (X) ] ^  <<  I ^ [ Li 2 (X) ] ^  for  the  case  of  nearly  complete 

modulation  and  small  ground  state  perturbations  since  [ Li^ ] j/ [ Li2 ] Q  ~ 

10  Also  the  term  [Li2(B) ] q [Li ]  ^  <<  [ Li2 (B)  ]  ^  [ L.i ] q  since  [ Li2 (B)  ]  0 

=  [ Li2 (B) ] !  for  nearly  complete  modulation  and  [Li]/[Li]0  <•  10  ^  based 

on  the  previous  steady  state  analysis  which  showed  that  [Li*]/ [Li]  £ 

-4 

10  .  By  eliminating  these  negligibly  small  terms  it  is  seen  that: 

(iu)  +  A2  +  K2  [  Li  ]  o  +  K6[Ar]0)[Li2(B)]1e1<*lB  =  IjCJ/hvtLi^X)  ]o 
which  reduces  to 

tan  <)>B  =  uj/(A2  +  A^  +  K2  [Li ]q  +  K6[Ar]Q) 

by  use  of  the  Euler  identity. 

Although  this  analysis  is  somewhat  more  complicated  than  for 
collisionally  unperturbed  emission,  the  results  are  nonetheless  of  the 


T 


same  form;  tan  <J)  =  U)T  where  1/T  =  [ Li  1 Q  +  [ Air ] q  .  Here  T 

is  easily  recognized  as  the  kinetically  controlled  lifetime  of  Li^CB) 
due  to  both  radiative  decay  and  collisional  losses  of  the  excited 
state.  The  results  of  the  phase  shift  experiments  can  therefore  be 
analyzed  in  the  standard  Stern-Vollmer  form  where  1/T  or  to/tan  ($>  is 
plotted  versus  pressure  or  species  densities.  The  Stern-Vollmer 
intercept  is  then  the  radiative  decay  rate  in  the  absence  of 
collisions  and  the  slope  of  the  line  is  the  rate  constant  for 
collisional  loss  of  the  excited  state  species. 

The  derivations  of  the  phase  shift  expressions  for  the  other 
excited  species  present  in  the  Li  only  and  Li/Mg  systems  are  presented 
in  Appendix  F  and  summarized  in  Figure  30.  For  simplicity  of 
presentation  the  "kinetic  decay  rates"  A^ 
in  the  equations  for  the  phase  shifts  and  are  separately  defined.  In 
all  of  the  derivations  the  same  previously  defined  assumptions  are 

2 

used.  The  expressions  for  the  phase  shifts  between  L^CB)  and  Li(  P) 
as  well  as  L^CB)  and  LiMg(A)  yield  the  same  form  as  above.  The 
expression  for  the  phase  of  L^fB)  relative  to  L^CA)  is  much  more 


through  A^,.  have  been  used 


complicated  due  to  the  existence  of  multiple  formation  processes 
involving  several  excited  state  species  each  having  a  unique  phase 
relationship  with  the  laser  (i.e.  second  order  kinetics). 


— — MnRPPPiiwMi  a  1  ii  1 1  i^B 

p  r  ^ 


I.  Li2  (B)  -  Li(2P)  Phase  Shift  No  Mg: 
u)  /tan^j  =  AK1 

AK]  £  Aj  +  lLi2(X)]K3  +  [Ar]K9 

II.  Li2(B)  -  Li2(A)  Phase  Shift  No  Mg: 

w  /tan<}>3  =  (AK3A4  +  (AK3AKl+U)2>A)/(A4  +  ^AK1+AK3^ 

AK3  h  A3  +  lAr]K? 

III.  Li2 (B)  -  LiMg(A)  Phase  Shift: 
a)  /tan05  =  AK5 

A]<5  =  A5  +  (Li]K4  +  [Ar)(K5  +Kg) 


Figure  30  -  Li/LiMg  Phase  Shift  Expressions 


i 
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C.  Experimental  Configuration 


1.  '.eady  State  Efficiency  Measurements 

Measurements  were  made  of  the  effect  on  emission  of  addition  of 
Mg  to  a  heatpipe  containing  just  Li  when  excited  by  an  Ar  ion  laser. 
Changes  in  emission  intensity  from  the  Li^CB-X)  and  Li^CA-X) 
transitions  were  interpreted  as  resulting  from  kinetic  processes 
involving  Mg.  Additionally,  by  comparing  the  change  in  Li^CB) 
emission  upon  Mg  addition  with  the  appearance  of  LiMg*  emission  the 
efficiency  of  excimer  production  was  determined. 

The  experiments  were  performed  using  the  same  configuration  as 
for  the  initial  laser  chemiluminescence  studies  described  in  Chapter  2. 
Emission  spectra  were  recorded  from  450nm  to  llOOnm  with  only  Li 
present  and  with  both  Li  and  Mg  present.  A  bandpass  filter  was  used 
to  exclude  the  second  order  Li^CB)  emission  which  obscured  the  LiMg 
emission  in  the  region  of  940nm- lOOOnm .  The  areas  under  the  emission 
curves  were  then  integrated  after  correcting  for  variations  in  the 
detector  sensitivity  over  the  wavelength  range,  as  well  as  the 
bandpass  filter  response.  The  relative  intensities  and  changes  in  the 
intensities  of  both  LiMg*  and  Li^CB)  were  then  used  with  the  steady 
state  equations  from  Figure  28  to  relate  various  kinetic  rates. 

2.  Phase  Shift  Measurements 

Time  delays  (or  phase  shifts)  of  one  optical  signal  relative  to 
another  were  recorded  at  various  frequencies  and  various  species 
number  densities(temperatures) .  As  an  example,  the  phase  delay  of 
the  LiMg  emission  relative  to  the  L^CB)  emission  was  measured  and 
related  via  the  phase  shift  equations  shown  in  Figure  30  to  the  LiMg 
decay  rate. 
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The  phase  shift  experiments  were  performed  using  a  slight 
modification  of  the  laser  chemiluminescence  experimental  configuration 
shown  in  Figure  11.  The  mechanical  chopper  was  replaced  by  an 
electro-optic  modulator  which  could  be  operated  at  frequencies  between 
lOKHz  and  1MHz.  Additionally,  a  high  frequency  lockin  amplifier  was 
required  to  extend  the  operating  range  to  1MHz.  Since  according  to 
tan<J>  =  U)T  it  is  expected  that  a  phase  shift  of  .35°  will  be  obtained 
at  a  modulation  frequency  of  lOOKHz  for  T  =  lOnsec,  it  is  necessary  to 
use  higher  modulation  frequencies  to  study  short  lived  species.  In 
this  series  of  experiments  the  phase  relationships  between  the  laser, 
Li^CB),  Li*,  Li^CA)  and  LiMg  were  measured  both  with  and  without  Mg 
present  at  various  species  densities. 

3.  Intensity  Measurements 

An  additional  set  of  measurements  were  made  where  the  intensity 
of  L^CB)  fluorescence  was  recorded  versus  temperature  in  the  absence 
of  Mg.  Since  radiation  trapping  of  Li2(B)  is  negligible  and  the 
absorption  of  the  laser  linearly  tracks  the  L^CX)  concentration  any 
deviation  from  proportionality  between  fluorescence  intensity  versus 
temperature  and  the  vapor  pressure  curves  must  represent  collisional 
deactivation  of  Li2(B)  by  Li  atoms.  The  experiment  was  performed  in 
the  same  configuration  as  shown  in  Figure  11  where  instead  of  covering 
the  entire  spectral  range,  only  the  Li^CB)  intensity  as  a  function  of 
temperature  was  monitored. 
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D.  Results  and  Analysis 
1 .  Introduction 

The  previously  developed  steady  state  rate  equations  were  used  to 
interpret  the  experimentally  determined  steady  state  efficiency,  phase 
shift  and  intensity  data.  The  results  were  interpreted  for  two 
separate  cases:  Li  alone  present  and  Li  present  in  combination  with 
Mg.  The  Li  only  case  is  analyzed  first,  since  the  rate  constants 
determined  for  this  case  provide  essential  inputs  to  the 
interpretation  of  the  combined  Li/Mg  system.  Because  of  the 
complexity  of  the  analysis  extensive  use  is  made  of  appendices  for 
discussion  of  ancillary  issues. 

2.  Li  Only  System 

a.  Determination  of  -  Li^CX)  +  Li(2P)  ■*  Li^(A)  +  Li(^S) 

Li (2S) 

Figure  31  shows  a  Stern  Vollmer  plot  derived  from  the  Li^CB)  - 

2 

Li (  P)  phase  shift  data  where  the  data  was  related  to  the  equation 


u)/tan<j>  =  AT  +  [Ar]K9  +  [ Li 2 (X)  ] K3 


Each  data  point  is  an  average  of  nine  individual  points;  three  each  at 

lOOKHz,  150KHz,  and  200KHz.  The  error  bars  show  one  standard 

deviation  of  the  data.  Since  during  the  course  of  the  experiments  the 

Ar  density  was  held  constant  and  only  the  Li^CX)  density  was  varied, 

the  slope  of  the  curve  can  only  be  attributed  to  rate  constant  K^,  the 

2 

rate  of  production  of  Li^CA)  from  collisions  of  Li^CX)  with  Li (  P), 

- 10  3 

equal  6.5  x  10  cm  /s.  The  intercept  is  the  sum  of  the  remaining  two 


constant  terms.  The  intercept  was,  therefore,  used  to  establish  an 

-13  3 

upper  limit  to  the  rate  constant  <  5  x  10  cm  /s.  Correspondingly 
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if  the  quenching  term  is  considered  negligible  then  an  upper  bound  of 
2 p  sec  is  obtained  for  the  radiation  trapped  lifetime.  This  is  in 

reasonable  agreement  with  the  value  of  10(J  sec  obtained  for  the  Na 

„  48 

system. 

b.  Determination  of  K2  -  Li2(B)  +  Li(2S)  -»  Li2(X)  +  Li(2P) 

A  determination  of  K2  was  obtained  from  measurements  of 
the  phase  shift  between  Li2(B)  and  Li2(A)  at  85Q°C  and  895°C. 

The  data  is  shown  in  Table  18. 


T(°C) 


ui/tan  0 


3.11.9  X  10 


6.81.3  X  10 


Table  18  -  Li2(B)  -  Li2(A)  Phase  Shift  (No  Mg)  Phase  Shift 


Again  each  entry  for  w/tan0  was  an  average  of  9  values;  3  each 
taken  at  lOOKHz,  150KHz,  and  200KHz.  The  uncertainties  shown  are 
the  standard  deviations  of  the  data.  The  data  were  related  to  the 
expression  shown  below  as  described  in  Appendix  F,  section  II. 


tanOyw  =  (A4  +  (AK1  l  AK3)A)/(AK3A4  +  (A^  +  u,  )A) 


where 


A  =  [Li2(X)]0K3(Li]0K2/(uj2  +  A2j) 


Rate  constant  K2  was  determined  to  be  2.711.3  x  10  cm3/s  over  the 
range  of  850°C  -  895°C. 

A  second  determination  of  K2  was  made  from  the  results  of  the 
experiments  where  Li2(B)  fluorescence  intensity  was  recorded  as  a 
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function  of  temperature  as  shown  in  Figure  32.  The  data  is  related  to 
the  expre'sion  shown  below  as  derived  from  the  steady  rate  equations 
in  Appendix  H: 

y  =  1  +  K2  [Li]/A2 

where  y  is  the  ratio  of  the  recorded  Li2(B)  emission  intensities  at 
temperatures  and  times  the  ratio  of  [Li2(X)l  at  temperatures  Tq 
and  T.:  y  =  ( [Li2(B) lTi/ [Li^(B) ] To )x ( [Li2(X) )Tq/ t Li 2 (X ) ] T± ) •  In  the 
absence  of  collisions  therefore  y=l.  With  the  experimental  values  of 
y,  the  values  of  [Li]  and  [Li  (X)]  at  a  given  temperature  from 
Appendix  B,  and  the  Li2(B)  A  coefficient  from  Appendix  E,  an  average 
value  of  K2  over  the  temperature  range  of  the  data  was  determined  as 
shown  in  Appendix  H.  K2  (the  rate  of  deactivation  of  Li2(B)  by  Li) 
was  determined  to  be  4.3±1.5  x  10  ^  cm3/s  over  the  temperature  range 
of  800°C  -  9,  JC. 

The  values  of  rate  constant  K2  obtained  from  two  different 
experiments  were  within  experimental  uncertainties.  However  both 
analyses  required  the  use  of  estimates  for  unknown  physical 
parameters.  The  intensity  data  analysis  used  the  estimated  value  of 
the  Li2(B)  A  coefficient  from  Appendix  E.  The  phase  shift  analysis 
required  rate  constant  which  was  previously  determined  as  well  as 
the  estimated  Li2(B-A)  radiative  rate  also  from  Appendix  E.  The  value 
of  K2  used  in  future  analysis  will  be  the  intermediate  value  of  the 
two  measurements:  K2  =  3.5  ±2.1  X  10  ^cm3/s. 

The  rate  constants  determined  for  the  Li  only  system  are 
summarized  in  Table  19.  Utilizing  the  values  of  rate  constants  K2  and 
expressions  can  be  developed  for  the  kinetic  lifetimes  of  Li(2P) 
and  Li2(B),  A^  and  A^j  respectively  by  using  the  algebraic 


expressions  as  given  in  Figure  25.  The  rates  that  were  determined  for 


the  Li/Li^  system  were  of  the  same  order  of  magnitude  as  those  for  the 
Na/Na^  system.  This  result  which  was  expected  lends  confidence  to  the 
experimental  approach  that  has  been  employed  to  determine  the  kinetic 
rates. 


Table  19  -  Li  System  Kinetic  Rates 

K2  =  3.5  +  2.1  X  10_1°  (cm3/s) 

_  10 

K3  =  6.5  ±  2.5  X  10  cm3/s) 

.13 

K9  <  5  X  10  cm3/s) 

A1  (Radiation  Trapped)  <  2jJ  sec 

Ar1  =  4.5  X  105  +  [ Li 2 (X) ]  6.5  X  10_1°  (1/s) 

A^  =  1.0  X  108  +  l Li ]  3.5  X  10_1°  (1/s) 

3.  Li/Mg  System 

The  majority  of  the  Li/Mg  system  rate  constants  were  derived  from 
measurements  of  the  efficiency  of  light  production  experiments.  A 
summary  of  these  results  are  shown  in  Table  20.  All  of  the  efficiency 
measurements  were  made  at  a  temperature  of  890°C  with  100  torr  of  Ar 
present  using  one  watt  of  488  nm  Ar  ion  laser  light  for  excitation. 

The  integrated  emission  intensities  of  L^tB),  L^tA)  and  LiMg(A)  were 
determined  both  with  and  without  Mg  present.  The  recorded  intensities 
were  corrected  for  wavelength  dependent  detector  and  filter  responses. 
The  ratios  of  the  integrated  intensities  of  the  various  emissions  are 


recorded  in  Table  20. 


Table  20  -  Intensity  Ratios 


T  =  890°C,  P(Ar)  =  lOOtorr,  A  =  488nm 


I.  Li  Only 

a.  Li2(A)/Li2(B)  =  .52 

II.  Li  +  Mg 

a.  Li2(A)/Li2(B)  =  .63 

b.  LiMg/Li2(B)  =  .09 

e.  LiMg/Li^CA)  =  .14 

III.  Changes  in  Intensity 

a.  Li^CB)  with  Mg/L^CB)  without  Mg  =  .35 

b.  Li2(A)  with  Mg/Li2(A)  without  Mg  =  .42 

IV.  LiMg  Conversion  Efficiency 

LiMg/(Li2(B)  without  Mg  -  Li2(B)  with  Mg)  =  .048 

a.  Determination  of  +  K^q  -  Li2(B)  +  Mg  ■*  LiMg(A)  +  Li 

Li2(A)  +  Mg 

The  total  rate  of  deactiviation  of  Li2(B)  by  collision  with  Mg 
atoms  was  determined  from  the  data  in  Table  20.  According  to  the 
kinetics  model  Mg  may  either  quench  Li2(B)  to  Li2(A)  (Kjq)  or  react  to 
form  LiMg(A)  (Kj).  Experimentally,  a  65%  reduction  of  Li2 (B)  emission 
was  measured  after  the  addition  of  Mg,  which  is  indicative  of  the  fact 
that  the  Mg  deactivation  by  collision  is  two  times  faster  than 
radiative  decay.  This  experimental  result  can  be  related  to  the 
steady  state  equations  for  Li2(B)  with  and  without  Mg  present.  The 
ratio  of  the  two  steady  state  equations  taken  from  Figure  21  is  shown 
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below.  The  subscripts  P  and  M  refer  to  the  pure  Li  case  and  the  mixed 
Li  +  Mg  case  respectively. 


[Li2(B)]p/[Li2(B)]M  =  (Aj£2  +  (Kj  +  K10)Mg)/AK2 

By  using  the  expression  for  A ^  from  Table  19  together  with  the  values 
of  [Li^]  =  7.9  x  lO^/cm^,  [Mg]  =  8.1  x  10^/cm^  from  Appendix  B  and 

[Li2(B) ] p/ [Li2(B) =  2.86  from  Table  20,  then  (K^  +  ^jo)  e9ua^  2.9  -  1 

lrt-10  3 
x  10  cm  . 

b.  Determination  of  -  LiMg(A)  +  Li  Li2(A,a)  +  Mg 
In  addition  to  measurements  of  the  light  conversion  efficiency, 
the  phase  shift  between  LiMg(A)  and  Li  (B)  was  measured  for 
temperatures  of  880°C  -  910°C.  The  results  are  shown  in  Table  21. 

Each  measurement  is  the  mean  of  at  least  seven  individual 

Table  21  -  LiMg  Phase  Shift  Data 
A  =  488  nm,  P(Ar)  =  100  torr 


Temperature 

u)/tan({» 

880°C 

+ 

3°C 

1.4 

+ 

.5 

X 

10 

900°C 

+ 

3°C 

.5 

+ 

.2 

X 

10 

905°C 

± 

3°C 

.6 

+ 

.3 

X 

10 

910°C 

± 

3°C 

1.5 

+ 

.6 

X 

10 

900°C 

± 

t— * 

o 

o 

n 

1.0 

+ 

.5 

X 

10 

measurements  spanning  the  frequency  of  50  -  250  KHz.  Because  the 
signal  to  noise  ratio  was  not  high,  the  data  tends  to  have  rather  large 
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i 


standard  deviations. 


The  LiMg  phase  shift  data  can  be  related  to  the  phase  shift 
equation  from  Appendix  F: 

AK5  =  u)/tan<J>5  =  A5  +  [Li  ]K^  +  [Ar](Kg  +  K^) 

where  A^  is  the  LiMg*  A  coefficient,  [Li]  is  the  deactiviation  rate 

from  Li  collisions  with  LiMg(A)  and  (Kg  +  K^)  are  Ar  deactivation  rate 

constants.  The  sum  of  rate  constants  Kg  +  is  determined  in 

-12  3 

Appendix  H  to  be  <  2  x  10  cm  /s  based  on  data  from  Chapter  2, 

Figure  19,  where  the  LiMg  emission  intensity  was  measured  as  a 

17  3 

function  of  Ar  pressure.  For  100  torr  (8.3  x  10  /cm  )  of  Ar  the  term 

[Ar](K0  +  Kc)  <  1.6  x  10^  which  is  negligible  in  comparison  to  the 
o  5 

measured  value  for  AR5  of  1  x  10  .  Therefore,  the  sum  A^  +  [Li)K^  is 

approximately  1  x  10^/s.  Since  theoretically  A^  'v  1.3  x  10^/sec  as 

shown  in  Appendix  E,  its  contribution  to  A^  is  also  negligible. 

Therefore,  A^^  reduces  to  a  single  term  [Li]K^  =  10^/s,  and  thus  = 
-10  3 

1.1  +  .5  x  10  cm/s.  The  last  approximation  will  be  further 
reinforced  by  showing  that  experimentally  A,.  •v  4  x  10^/sec. 

c.  Determination  of  -  L^CB)  +  Mg  •*  L^CA)  +  Mg 
The  efficiency  experiment  data  in  Table  20  shows  that  Mg  inter¬ 
acts  with  L^CB)  in  such  a  way  as  to  produce  an  enhanced  density  of 
Li2(A),  since  the  L^CAVL^CB)  ratio  is  .52  without  Mg  and  .63  with 
Mg.  The  kinetic  model  accounts  for  this  result  in  two  ways.  Mg  may 
directly  quench  Li2(B)  to  Li2(A)  -  (K1Q)  or  Mg  may  react  with  Li2 (B) 
producing  the  excimer  LiMg(A)  -  (Kj)  and  then  the  excimer  molecules 
may  in  turn  collide  with  Li  atoms  producing  L^CA)  -  (K^). 


An  upper  bound  to  K^q  can  be  obtained  by  assuming  K^a=0  and  by 
then  comparing  the  Li^(A)  steady  state  equations  with  and  without  the 
Mg  terms.  The  resulting  ratio  is  shown  below  where  the  pure  Li  case 
is  subscripted  P  and  the  mixed  Li/Mg  case  is  subscripted  M: 

(Li2(A)]M  [Li(2P)]M(Li2(X)]K3  +  [Li2(B)]M(A4  +  \WlQ) 

[ Li 2 ( A ) ]p  [Li(2P)]p[Li2(X)]K3  +  [Li2(B)]pA4 

Appendix  I  develops  the  relationship  [ Li (  P)]  =  10[Li2(B)].  Making 

this  substitution  above  and  using  the  ratios  (Li2(A) 1^/ (Li2(A) ]p  = 

.42,  [ Li2 (B) ] M/  [Li2(B)]p  =  .35  from  Table  20,  [ Li2 (X) ]  and  [Mg]  equal 

4.1xl0^^/cm^  and  8.1xl0^^/cm^  at  890  C  respectively,  A^  equal  lxlO7 

-10  3 

from  Appendix  E  and  K3  equal  6.5x10  cm  /s  from  section  2. a,  then  J<10 

<  9.1x10  12  cm'Vs.  Thus,  the  value  of  is  less  than  3%  of  the  sum 

of  as  determined  in  section  3. a.  Therefore,  the  calculation 

in  this  section  substantiates  that  is  the  rate  constant  for  the 

dominant  Li2(B)  deactivation  process  and  is  equal  to  2.9  x  10  ^ 

3, 

cm  /s . 

d.  Determination  of  LiMg(A)  A  Coefficient  -  A,- 
The  value  of  A^ ,  the  LiMg(A)  A  coefficient  can  be  evaluated  from 
a  combination  of  the  LiMg  phase  shift  data  and  relative  intensity 
data.  The  steady  state  rate  equation  for  LiMg  (Figure  21)  can  be 
rewritten  as 

[LiMg(A)  ] /  ( Li 2 (B)  ]  -  [MgjK^A^ 


The  ratio  on  the  left  hand  side  of  the  equation  is  obtained  from  the 


product  of  the  integrated  intensity  ratio  I (LiMg(A) )/I (Li2(B) )  =  .09 
from  Tab-  20  and  the  ratio  of  the  Li2(B)  and  LiMg(A)  A  coefficents 

g 

(10  /A^).  The  right  hand  side  can  be  evaluated  since  A^  was 

previously  determined  to  be  1  x  102/s  from  the  phase  shift  data,  (Mg] 

at  890°C  is  8.1  x  10^/cm^  from  Appendix  B.  Now  since  >>  (i.e. 

that  the  chemical  reaction  rate  of  Mg  with  Li^(B)  is  much  larger  than 

-10  3 

the  quenching  rate)  and  Kj  =  2.9  +  1  x  10  cm  /s,  the  resulting 

value  of  A^  =  3.8  x  10^/s.  The  radiative  lifetime  for  LiMg(A)  is 

therefore  2 . 6(j  seconds.  The  experimentally  obtained  value  of  the  LiMg 

A  coefficient  is  between  the  two  theoretical  values  1.3  x  10^/s  for 

4 

the  bound-free  emission  and  .2  x  10  /s  for  the  bound-bound  emission  as 
shown  in  Appendix  E. 

e.  Determination  of  K^a  (Formation  of  Li2(A)  from  LiMg) 
and  (Formation  of  Li2(a)  from  LiMg) 

If  now  it  is  assumed  that  all  the  additional  Li^CA)  emission  that 
occurs  upon  addition  of  Mg  is  due  to  reaction  of  Li  with  LiMg(A)  (K^a) 
and  the  formation  of  Li2(A)  due  to  the  quenching  reaction  is 
negligible  (K^q=0) ,  then  an  upper  bound  on  K^a  can  be  determined  from 
the  same  Li2(A)  data  used  above.  As  shown  above: 


lLi2(A)]M 

[Li2(A)]p 


[Li(2P)]M(Li2(X)]K3  +  (LiMg(A) ]  [Li ]K4a  +  (Li^B)]^ 
[Li(2P)]p[Li2(X)]K3  +  [Li2(B)JpA4 


and  using  the  same  values  as  before  plus  the  ratio  (LiMg(A) ] / (Li2(B) ] ^ 

=  (.09)/ (.35)  x  (A2/Ar)  =67.7  from  Table  20  and  Appendix  E,  then  K^a 
- 1 3  3 

<  4.8  X  10  cm  /s.  The  value  of  K,  is  thus  found  to  be  small 

4a 

compared  to  the  total  value  calculated  for  K,  in  section  3.b. 


-10  3 

Therefore,  =  1.1+.5  x  10  cm  /s  implies  that  the  reaction  of 

the  excimers  with  Li  atoms  yield  predominantly  Li^fa)  rather  than 
Li2(A). 

Finally,  since  the  production  of  Li2(A)  may  result  from  either 

direct  quenching  from  the  B  state  or  reactions  with  LiMg  and  the 

balance  of  the  two  cannot  be  determined  from  the  data,  only  upper 

bounds  for  the  two  rate  constants  K,„  and  K,  can  be  determined. 

10  4a 

A  summary  of  the  rates  determined  for  the  Li/Mg  system  is  shown 
below: 

K1  =  2.9  ±  1  x  10'10cm3/s 

K,  <  4.8  x  10  ^3  cm3/s 
4a 

-10  1 

K4fa  =  1.1  ±  .5  x  10  cm  /  s 

A5  =  3.8  x  105/s 

K1q  <  9.1  x  10-12  cm3/s 

Kg  +  K5  <  1.3  x  10'12  cm3/ s 

Figure  26  summarizes  the  specification  of  the  LiMg  model.  Either  a 
measurement  has  been  made  or  an  upper  bound  determined  for  every  rate 
in  the  model. 


107 


Table  22  -  LiMg  Model  Rate  Measurements 


A.  Excitation 


Li2(X)+hi>-»Li2(B) 
B.  Primary 


Li2(B)+Mg-»LiMg(A)+Li 

K1(cm3/s) 

2.9+lxlO*10 

Li2(B)+Li-»Li2(X)+Li(2P) 

K2(cm3/s) 

3.5±2.5xl0*10 

Secondary 

Li(2P)+Li2(X)^Li2(A,a)+Li 

K3(cm3/s) 

6.5±2.5xl0*10 

LiMg(A)+Li->Li2(A,a)+Mg 

K4(cm3/s) 

<4.8xl0*13;  l.l+.5xl0*10 

LiMg(A)+Ar-»Li(2P)+Mg+Ar 

K5(cm3/s) 

-12 

<2x10 

Radiation 

Li(2P)->Li+hv 

Aj  (1/s) 

3 . 72xl07 

Li2(B)-*Li2(X)+hv 

A2(1/s) 

g 

1x10  (Appendix 

Li2(A)->Li2(X)+hv 

A3(1/s) 

5.56xl07 

Li2(B)^Li2(A)+hv 

A4(1/s) 

lxl07(Appendix 

LiMg(A)->LiMg+hv 

A5(1/s) 

3.8xl05 

Quenching 

Li2(B)+Ar-*Li2(X)+Ar 

K6 

<3x10* 13 

Li2(A)+Ar-+Li2(X)+Ar 

K7 

<10  ^(Assumed) 

LiMg(A)+Ar->LiMg(X)+Ar 

K8 

<2xl0*12 

Li(2P)+Ar  ,Li-»Li+Ar  ,Li 

K9 

<5xl0*13 

Li2  (B)+Mg-»Li2  (A)+Mg 

K10 

- 12 

<9.1x10 
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Chapter  IV  -  Summary,  Conclusions  and  Recommendations 


During  the  thesis  work,  two  new  diatomic  species  (LiMg  and  LiCa)  were 
observed  and  characterized.  Using  laser  induced  chemiluminescence,  a 
characteristic  emission  spectra  was  obtained.  Ab  initio  calculations 
were  used  to  construct  a  synthetic  spectrum  which  was  in  good  agreement 
with  the  experimental  data.  Additionally,  the  formation  and  decay 
kinetics  for  the  LiMg  system  were  examined  theoretically  and  experi¬ 
mentally.  It  was  concluded  that  LiMg  was  formed  in  an  excited  state 
as  a  result  of  reaction  of  excited  L^CB)  with  Mg  atoms.  The  principle 
decay  channels  included  reaction  back  to  Li^CA)  by  collisions  of  the 
LiMg  excimer  with  Li  atoms  and  radiative  decay. 

Several  additional  experiments  could  be  performed  which  would  add 

to  the  body  of  knowledge  of  the  Group  IA/IIA  excimers.  Other  combinations 

of  IA/IIA  metals  could  be  studied  in  chemiluminescence  experiments  using 

experimental  apparatus  in  which  the  two  metal  pressures  could  be  varied 

independent  of  one  another.  A  vapor  pressure  mismatch  between  Na  and 

Ca  appeared  to  be  the  major  reason  for  not  observing  emission  from  the 

NaCa  excimer.  Additionally,  the  kinetics  of  the  LiCa  system  could  be 

studied.  Finally,  further  experiments  could  be  performed  on  the  LiMg 

3 

system.  The  reaction  of  Mg(  P)  with  Li^CX)  resulting  in  LiMg(B)  +  Li 
could  be  studied  and  the  data  used  to  better  characterize  the  LiMg 
potential  energy  curves.  Also  the  potential  of  LiMg  as  an  optically 
pumped  laser  candidate  could  be  examined. 
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Appendix  A  -  Experimental  Equipment 


Equipment 

GaAs  Photomultiplier  Tube 
Housing  &  Cooler 
Ge  Detector 

Scanning  Monochromator 

Low  Frequency  Lockin  Amplifier 

High  Frequency  Lockin  Amplifier 

Preamplifier 

Chart  Recorder 

Chopper 

Electro-optic  Modulator 

Function  Generator 

Universal  Counter 

White  Light  Source 

Ar  Ion  Laser 

CW  Dye  Laser 

Oscilloscope 

Power  Meter 

Pulsed  Dye  Laser 

High  Frequency  Preamplifier 


Manufacturer  &  Model 
RCA  C31004A 

Products  for  Research  TE1777S 
Applied  Detector  Corp.  //304 
McPherson  2051 

Princeton  Applied  Research  124A 
Princeton  Applied  Research  5202 
Princeton  Applied  Research  181 
Hewlett  Packard  9313 
Princeton  Applied  Research  125A 
Coherent  Associates  317 
Tektronix  FG  504 
Tektronix  DC  505A 
Unknown 

Coherent  Radiation  CR-8 
Coherent  Radiation  590 
Hewlett  Packard  1744A 
Scientech  362 
Phase-R  DL-2100B 
Hewlett  Packard  461A 
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Appendix  B 


The  species  number  densities  used  in  the  kinetics  calculations  were 
determined  utilizing  data  from  Table  II  in  Nesmeyanov^  which  supplied 


the  fitting  coefficients  for  L^,  Li,  Mg  and  Ca  to  the  equation  log 

P  =  A  -  B/T  +  CT  +  D  log  T  as  shown  below.  The  calculated  pressures 

3 

were  then  converted  to  number  densities  (1/cm  )  assuming  ideal  gas 
behaviour . 


Li 

Li2 

Mg 

Ca 

A 

10.34540 

18.37849 

16.79743 

53.65151 

B 

8345.574 

11139.62 

7844.230 

12826.75 

C 

-.0000884 

.0001634 

.0002548 

.0020016 

D 

-0.681060 

-3.032090 

-2.727980 

-14.48467 

Temp (°C) 

Li 

Li2 

Mg 

Ca 

650 

1.67E15 

3.10E13 

2.92E16 

4.72E14 

700 

4.42E15 

1.06E14 

6.76E16 

1.36E15 

750 

1.06E16 

3.21E14 

1.43E17 

3.47E15 

800 

2.31E16 

8.68E14 

2.81E17 

8.01E15 

850 

4.71E16 

2.13E15 

5.16E17 

1.70E16 

870 

6.13E16 

2.99E15 

6.47E17 

2.24E16 

875 

6.55E16 

3.24E15 

6.84E17 

2.40E16 

880 

6.98E16 

3.51E15 

7.23E17 

2.57E16 

885 

7.44E16 

3.81E15 

7.63E17 

2.74E16 

890 

7.92E16 

4.13E15 

8.06E17 

2.93E16 

895 

8.44E16 

4.46E15 

8.50E17 

3.13E16 

900 

8.97E16 

4.83E15 

8.96E17 

3.34E16 

905 

9.54E16 

5.22E15 

9.44E17 

3.56E16 

910 

1.01E17 

5.63E15 

9.94E17 

3.79E16 

915 

1.08E17 

6.08E15 

1.05E18 

4.03E16 

920 

1.14E17 

6.55E15 

1.10E18 

4.29E16 

115 


Appendix  C  -  Determination  of  Mg£  Density  in  a  Pure  Mg  Gas 

The  pressure  for  a  nonideal  gas  can  be  expanded  in  a  power  series 

in  the  number  density: 

P  =  nkT[l  +  nB(T)  +  n2C(T)  +....] 

where  B(T)  and  C(T)  are  the  second  and  third  virial  coefficients.^1  The 

terms  in  the  power  series  represent  the  contribution  to  the  pressure  of 

the  monomers,  dimers,  trimers,  etc.  The  term  of  interest  for  Mg^  is 
2 

n  kTB(T).  B(T)  can  be  calculated  from  the  expression: 

B(T)  =  2rtJoR2dR(l-e‘U(-R)/kT) 

where  R  is  the  internuclear  separation  and  U(R)  is  the  potential  energy 

62 

curve  for  the  ground  state.  For  Mg^  the  U(r)  can  be  approximated  by  the 
square  well  shown  below: ^ 

U(r)  =  *  for  r  <  r  ,  r  =  3& 
o’  o 

=  -e  for  r  <  r  <  ar  ,  -£  =  300/cm 
o  o 

=  0  for  r  <  ar  :  ar  =  6$ 
o  ’  o 


then 

B(T)  =  2rtJ^o  R2dR(l-eC/kT)  +  2rtro3/3 

-  (27ir3/3)[  1  -  (a3  -  1  )(  ee/kT  -  1  )] 

-22  3 

for  a  temperature  of  900°C  then  B(T)  =  1.8  x  10  cm  and  using  the  value 
for  N(Mg)  at  900°C  then  N(Mg2)  =  1.5  x  10^/cm3  or  Mg2/Mg  =  1.6  x  10  . 

By  analogy,  B(T)  for  LiMg  is  a  square  well  100/cm  deep  and  width 
2A-6A,  B(T)  =  6.9  x  10  23  cm3  and  N(LiMg)  =  5.6  x  1012/cm3  the  ratio 
LiMg/Mg  is  6  x  10  ^ . 
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Appendix  D  -  Mg^  and  LiMg  Density  Determination  in  a  Li/Mg  Gas  Mixture 

A  chemical  reaction  can  be  analyzed  by  means  of  the  Gibbs-Duhem 

relationship  to  determine  one  of  the  product  densities  given  that  the 

63 

reactant  densities  and  other  product  densities  are  known.  The  Gibbs- 
Duhem  relationships  for  the  reactions: 

Li ^  +  Mg  ^  LiMg  +  Li 
m82  +  Li  <-  LiM8  + 

are 


exp[-BE(LiMg)/kT][m(Li2)m(Mg)}1-5  B(LiMg)  Qy(Li2)  [LiMg] [Li] 
exp[-BE(Li2)/kT][m(LiMg)m(Li)]1‘5  B(Li2)  Qv(LiMg)  [Li2J [Mg] 


and 

exp[-BE(LiMg)/kT][m(Mg2)  m(Li)]1'5  B(LiMg)  Qv(Mg2)  [LiMg]  [Mg] 
exp[-BE(Mg2)/kT][m(LiMg)  m(Mg)]1’5  B(Mg2)  Qv(LiMg)  [Mg2J  [Li] 


where  m(  )  are  the  masses,  B(  )  are  the  diatomic  rotational  constants, 

BE(  )  are  the  bond  energies  in  1/cm,  and  Q  (  )  are  the  vibrational 
64 

partition  functions  for  the  respective  diatomic:  =  I  exp(-u>e(v+^)/kT) . 

The  vibrational  partition  functions  depend  on  the  vibrational  spacing  uj^  and 


m 

U) 

Be 

Qv 

BE. 


summed  over  the 

vibrational 

levels  v. 

The  values 

equations 

are  shown  in  the 

table  below: 

.65,66,67 

Li 

Mg 

Li 

LiMg 

Mg? 

7 

24.3 

14Z 

31.3 

4876 

351.3 

65.7 

51.1 

.67 

.20 

.09 

2.8 

2.7 

o.26 

8900 

135 

400 
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Then  [LiMg]  =  3.9  x  10'6lLi2][Mgl/[Li] ,  [ Mg2 1  =  5 .9 [LiMg] [Mg] / (Li ] 
and  at  900°C  using  the  data  from  Appendix  B  the  results  are  shown 
below: 

[LiMg]  =  1.9  x  10n/cm3;  [LiMg]/ [Mg]  =  2  x  10"7. 
[Mg2l  =  1.1  x  1013/cm3;  [Mg2]/[Mg]  =  1.2  x  IQ-5. 
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Appendix  E  -  Summary  of  Radiative  Lifetime  Data  and  Calculations 

I.  Li^CB)  State  Lifetime  Estimate 

There  are  no  literature  values  for  the  L^CB)  state  lifetime. 
Therefore  the  lifetime  will  be  inferred  from  the  relative  lifetimes  of 

68-72 

the  Na  A  and  B  states  reported  in  the  literature  '  as  shown  below: 

Dimer  Li«  Na«  K 

State  ill 

A1!  18 . 0ns±2  12.2ns±.3  . 

B  1Fl  (10ns)*  7.0ns±.4  11ns 

"Inferred  from  Na2(B)/Na2(A)  lifetime  ratio. 

The  L^CB)  lifetime  estimate  reflects  the  total  decay  rate  of  this 
state  which  radiates  to  both  the  A1!  and  X1!  states.  Since  both  the 
B-A  and  A-X  transitions  are  states  having  approximately  the  same  shape 
to  their  potential  energy  curves  and  the  same  Rg  it  is  assumed  that  the 
Franck  Condon  envelopes  are  approximately  the  same.^  The  relative  rate 
into  each  state  is  therefore  proportional  to  the  transition  frequencies 

cubed.  Therefore  the  L^CB-X)  radiative  rate  will  be  assumed  to  be 

8  7 

10  /s  and  the  L^CB-A)  rate  will  be  assumed  to  be  10  / s. 

An  additional  estimate  of  the  L^CB)  lifetime  was  obtained  by 

comparing  the  area  under  the  absorption  curves  from  Figure  12.  By 

comparing  the  B-X  and  A-X  absorbed  areas  a  ratio  of  the  (B-X)  to  (A-X) 

A  coefficients  was  obtained.  Using  the  literature  value  for  the 

L^CA-X)  A  coefficient  a  value  of  12nsec  ±  4nsec  was  obtained  for  the 

L^CB-X)  radiative  decay  rate.  The  large  uncertainty  in  this  estimate 

is  due  to  the  difficulty  of  separating  Li  atom  absorption  from  the 

overlapping  dimer  absorptions. 
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II.  Summary  of  LiMg(A2fi)  State  Lifetime  Calculations 


Estimates  of  the  bound  free  and  bound  bound  transition 
probabilities  were  calculated  by  Michels  using  the  potential  energy 
curves  shown  in  Figure  4  and  dipole  moment  functions  from  Figure  6. 
Since  the  dipole  moment  function  is  changing  rapidly  with 
internuclear  separation,  the  A  coefficient  also  changes  with 
vibrational  level.  The  majority  of  the  observed  emission  emanates 
from  the  v'  =  3,4,5  of  the  A  state.  The  average  A  coefficient  for 
these  three  levels  is  . 13  x  10^/s. 


Bound-Free  Bound-Bound 


v' 

A(l/s) 

v1 

v" 

A(l/s) 

0 

.6  x  lof 

0 

0,2 

.2,  1 

1 

.7  x  10* 

1 

0,2 

3,  10 

2 

.9  x  10° 

2 

0,2 

30,  100 

3 

.11  x  10^ 

3 

0,2 

200,  500 

4 

.13  x  lo' 

4 

0,2 

700,  2000 

5 

.15  x  10 ' 

5 

0,2 

3000,  7000 

6 

.17  x  10/ 

6 

0,2 

10000,  30000 

Private  Communication,  H.  Michels. 
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Appendix  F  -  Phase  Shift  Derivations 


I.  Li^tB)  -  Li(2P)  Phase  Shift  Derivation  (No  Mg) 

2 

The  phase  shift  of  Li(  P)  relative  to  Li^CB)  was  determined 

2 

using  the  rate  equation  for  Li(  P)  shown  below  from  Figure  27: 

d[Li*]/dt  =  [Li2(B)][Li]K2  -  [Li*](A1  +  [Li2(X)]K3  +  [Ar]Kg) 

Then  upon  Fourier  analyzing,  keeping  only  terms  of  frequency  w, 
ignoring  the  negligibly  small  cross  terms,  cancelling  the  e1U)t,  and 
algebraically  rearranging: 

[Li*]1exp(i<J»Li2(B)  -  i«J>Li*)  =  [Li2(B)  ]Q[Li] qK2 / 

(iiu  +  (Ax  +  [Li2(X))0K3  +  [Ar]QK9)) 

By  defining  A^  =  Ax  +  (Li2(X)]QK3  +  [Ar]0Kg 

♦j  =  <KLi2(B))  -  <KLi(2P)) 

A  e  [Li]0K2/(u»2  +  A22) 

And  then  substituting  into  the  equation  above: 

[Li  ]1exp(i4>1)  =  A[Li2(B)]1(iu)  +  AR1) 

Then  Re([Li*])  =  AlLi^B)]^ 

Im( [Li  ])  =  Atu[Li2(B)]1 
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Since  by  the  Euler  identity: 


[Li*] ^exp(i0j )  =  Re(|Li*])  + 


Thus : 


w/tan  4>1  =  AJ(1 


ilm( [Li*] ) 
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II.  Li 2 ( B )  -  Li  (A)  Phase  Shift  (No  Mg) 

The  phase  of  Li2(A)  relative  to  Li2(B)  was  determined  using 
the  rate  equation  for  Li^(A)  as  shown  below  from  Figure  27.  Since 
the  formation  of  Li^CA)  depends  on  both  the  direct  process  of  radiation 
from  the  B  state  and  a  second  order  process  where  Li*  is  formed  as 
an  intermediate  the  kinetics  become  quite  complex. 

d[Li2(A)]/dt  =  [Li"][Li2(X)]K3  +  [Li2(B)]A4-[Li2(A)](A3+K7[Ar]0) 

By  Fourier  analyzing,  keeping  only  the  terms  at  frequency  iu, 
ignoring  the  negligibly  small  terms  and  cancelling  the  exp(iwt) 
and  algebraically  rearranging: 

( Li2 (A) ]  1exp(i(j»(Li2(B) )-i<Ji(Li2(A)) )  =  ([Li  ]1exp(i0(Li*)-i0(Li2(B))) 
[Li2(X)]QK3  +  [Li2(B)]1A4)/(iu>  +  A3  +  K? [Ar]Q) 

Note  that  two  different  phase  shifts  remain  in  this  equation;  the 
phase  of  Li*  relative  to  Li2(B)  and  the  phase  of  Li2(A)  relative  to 
Li2(B).  By  making  the  following  substitutions: 

AK3  £  A3  +  K? [ Ar ] Q  ;  0(Li2(B)  -  Li2 (A) )  =  <J>3 
<t>j  5  0(Li*)  -  <|>(Li2(B)) 

03  H  0(Li2(B)  -  0(Li2(A)) 

And  rationalizing  the  denominator: 
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[  L  i  2 1 A )  ]  1exp(i4>3)  =  ([Li  J^expCi^)  [  L±2  (X)  )  QK3  +  [  Li  2  ( B  Jt  ]  jA^)  (iw+AK3)/ 


'"2  *  AK3> 


The  real  and  imaginary  parts  of  this  expression  are: 


Re([Li2(A)]1exp(i(J)3))  = 

(  [  Li2  (B)  ]  A^+Re  ( [  Li  ]  3exp(  itl)^^ ) )  [  Li2  (X) )  QK3)AK3+ujIm(  [  Li  ]  ^exp  ( if}^ ) )/ 
(u,2  +  AK32) 

Im([Li2(A)l  1exp(i([)3)  =  (uj(  [  Li  2  ( B )  ]A^  +  Re  ( (Li*)  ^exp^ )  [  L  i  2  ( X )  ]  ^ ) 

+  [Li2(X)]0K3AK3Im([Li*]1exp(i([)1))/(u)2  +  A^) 

Note  that  the  above  expressions  include  the  real  and  imaginary 
expressions  for  Li*  which  can  be  obtained  from  the  previous  section. 
Making  these  substitutions  and  applying  the  Euler  identity  results 
in  the  following  expression  for  4>3 . 

tan<l>3  =  w(A4  +  (AK1  +  AK3)A)/  (AK3A4  +  +  uj2)A) 


To  evaluate  the  equation  above  for  K2  using  the  measured  values  of 
U)/tan<()  requires  the  values  for  A^,  ^K1  ’  anc*  t  Li  2  (X)  ]  0^3  •  ^3 

assumed  to  be  just  A 3  since  a  typical  Ar  quenching  rate  of  10  ^ 
makes  the  quenching  term  in  A^  negligible.  A3  =  5.6  x  102/s 
from  Appendix  E,  A4  is  calculated  to  be  10^/s(see  Appendix  E),  the 
value  of  AV1  is  taken  from  the  phase  shift  data,  Table  19.  The 
numerical  values  for  the  different  variables  in  the  equation  are 
tabulated  below  along  with  the  data  for  u)/tan<}>  for  two  different 
temperatures : 


is 
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T(°C) 


A 


ui/tan0 


T(°C) 

2  2 
w  +aki 

[Li2(X)]K3 

aki 

A 

ui/tan0 

850 

4.1xl012 

1 .4xl06 

1 ,8xl06 

3 . 4xl0~7 

3.11.9x10' 

895 

1.2xl013 

2.9xl06 

3.4xl06 

2 . 4xl0_7 

6.81.3x10' 

By  using  the  appropriate  values  for  [ Li ] ^  from  Appendix  B,  is 

-in  -in  'l 

4.211.5x10  cm  /s  and  1.21.3x10  cm  /s  at  850°C  and  895°C  respectively, 

-10  3 

the  average  being  2.7x10  cm  /s.  Variations  in  temperature  do  not 
account  for  the  discrepancy  between  the  two  rate  constants. 

III.  LiMg  -  L^CB)  Phase  Shift  Analysis 

The  phase  of  LiMg  relative  to  L^CB)  was  determined  from  the 
rate  equation  for  LiMg  as  shown  below  from  Figure  27: 

d[LiMg(A)]/dt  =  -[LiMg(A)](fLi}K4+A5+[Ar](Kg+K5)  +  [Li2(B) ] [MgjKj 

Fourier  analyzing  ,  keeping  terms  only  at  the  fundamental 
frequency,  cancelling  negligible  terms  and  the  exp(iu>t),  and 
algebraically  rearranging  yields: 

C iui  +  [Li]QK4+A5+[Ar]  (Kg+K5))  (LiMg(A)  ]exp(i<t>^)  =  [Li ^ (B )  ] x  [Mg]QK1 

Defining  <)>,.  =  <}>(Li2(B)  -  LiMg(A)) 

AK5  =  [Li]K4  +  A5  +  [Ar] (Kg  +  K5) 

Substituting  into  the  equation  above  and  applying  the  Euler  identity: 

w/taruDs  = 
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Appendix  G  -  Li^CB)  Intensity  vs  Temperature 

Using  the  steady  state  solution  to  the  rate  equation  for  Li^CB) 
from  Figure  26,  an  expression  can  be  developed  for  the  temperature 
dependence  of  Li2(B)  emission  intensity.  Taking  the  ratio  of 
the  steady  state  Li^CB)  rate  equation  for  no  Mg  present  at  two  temper¬ 
atures  T.  and  T^  yields: 

1  f  3 

[Li2(B)]./[Li2(B)]f  =  ((A2  +  A4  +  K2[Li]f  +  [Q]K9)/[Li2(X)]f)/ 

((A2  +  A4  +  K2(Li],  +  [Q]Kg)/[Li2(X)].) 

If  the  initial  temperature  is  chosen  such  that  K2[Li]Q  =  A2+^4 

-11  A 

and  since  Kg  is  <<  10  cm  /s  then  for  100  torr  of  Ar,  +  A^  >>  (Ar]K 
thus  Kg  can  be  ignored.  Using  this  approximation  simplifies  the 
equation  above  to: 

Y  =  [Li2(B)]Q  [Li2(X))./[Li(B)].  [Li2(X)]Q  =  1  +  K2[Li]./A2  +  A4 
The  values  of  [Li2(X)  ]  ^./ [Li2  (X)  ]  ^  were  calculated  from  Appendix  B  and 

O 

data  points  from  844  C  to  893. 5°C  from  Fig  32  were  used  with  the 
equations  shown  above  to  determine  K2as  shown  below: 


T°(C)  X  K  (cm3/s) 


844 

1.15 

4.4 

X 

10-10 

849 

1.13 

3.3 

X 

10-10 

856.5 

1.17 

4.1 

X 

10-10 

864 

1.16 

3.4 

X 

10-10 

874 

1.23 

4.0 

X 

10-10 

883.5 

1.28 

4.2 

X 

10-10 

893.5 

1.40 

5.3 

X 

10-10 

-10  3 

The  value  of  determined  from  this  data  is  4.311.5  x  10  cm  /s. 
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Appendix  H  -  Ar  Quenching  of  LiMg(A) 


Using  the  data  from  Figure  18  in  Chapter  2  an  upper  bound  can  be 
obtained  for  the  sum  of  the  rates  of  deactivation  and  dissociation  of 


LiMg (A)  by  Ar(Kg 

+  k5). 

It  is  an  upper 

bound  since  the 

Li^CB)  was 

assumed  to  be  immune  to 

Ar  quenching. 

The  data 

is  summarized  below 

Total 

Ar 

Ar 

I 

l'lo 

K8  +  K5 

Pressure 

Pressure  Density 

425 

316 

2.6xl018 

132 

1.44 

1.7x10” 

344 

235 

1.9 

152 

1.26 

1.4 

280 

171 

1.4 

158 

1.2 

1.4 

204 

95 

7.8xl017 

170 

1.12 

1.5 

162 

53 

4.4 

185 

1.03 

.7 

122 

13 

1.1 

191 

1 

The  rate  constant  is  determined  by  taking  the  ratio  of  the  LiMg  rate 
equation  from  Figure  28  with  and  without  Ar  present: 

[LiMg(A)]A/[LiMg(A)]NA  =  (A5  +  [LijK4)/(A5  +  [Li JK^  +  [AR] (Kg  +  Kg) 

By  substituting  Ar  +  [LijK^  =  107/s  from  the  LiMg  phase  shift  data  from 
Table  21  and  I/Iq  as  the  ratio  of  the  LiMg  intensity  with  the  specified 
amount  of  Ar  to  the  LiMg  intensity  with  13  torr  of  Ar  (  Ar  effects 
assumed  negligible)  then: 

107((I/I0)  -  1)/ [Ar]  =  Kg  +  K5 

-12 

The  resulting  average  for  the  upper  bound  on  Kg  +  is  1 . 3  x  10 
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Appendix  I  -  [Li  ]  to  [Li^CB)]  Relationship 

An  expression  was  developed  for  the  concentration  of  Li*  in 
terms  of  Li^CB)  using  the  rate  constants  determined  in  section  2 
of  Chapter  3.  The  governing  steady  state  equation  from  Figure  28 
is  shown  below: 


[Li" ]  =  [Li2(B)][Li]K2/([Li2(X)]K3  +  Aj  +  [Ar]Kg) 

The  results  of  Chapter  3  section  2. a  showed  that  the  terms  A^  +  [Ar]K 
in  the  denominator  of  the  right  hand  side  of  the  equations  are 

O 

negligible  when  compared  to  [Li2(X)]K3.  Then  for  890  C 

[Li]  =  7.9  x  10^/cm^  and  [  Li2  (X)  ]  =  4.1  x  10^  from  Appendix  B,  and 

-10  3  -10  3 

K2  =  3.5  x  10  cm/s  and  =  6.5  x  10  cm  /s  from  Table  19, 


yields : 


Vita 
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